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Abstract 



It is well known that the s-wave Jost function for a potential, XV, is an entire 
function of A with an infinite number of zeros extending to infinity. For a repulsive 
V, and at zero energy, these zeros of the "coupling constant", A, will all be real and 
negative, A n (0) < 0. By rescaling A, such that A n < —1/4, and changing variables 
to s, with A = s(s — 1), it follows that as a function of s the Jost function has 
only zeros on the line s n = \ + ij n . Thus finding a repulsive V whose coupling 
constant spectrum coincides with the Riemann zeros will establish the Riemann 
hypothesis, but this will be a very difficult and unguided search. 

In this paper we make a significant enlargement of the class of potentials needed 
for a generalization of the above idea. We also make this new class amenable to 
construction via inverse scattering methods. We show that all one needs is a one 
parameter class of potentials, U(s; x), which are analytic in the strip, < Res < 1, 
Ims > T Q , and in addition have an asymptotic expansion in powers of [s(s — l)] -1 , 
i.e. U{s;x) = V (x) + gV 1 {x) + g 2 V 2 (x) + ... + 0(g N ), with g = [s(s - l)] -1 . The 
potentials V n (x) are real and summable. Under suitable conditions on the V^s and 
the 0(g N ) term we show that the condition, f °° \f (x)\ 2 V\(x)dx / 0, where f is 
the zero energy and g = Jost function for U, is sufficient to guarantee that the 
zeros g n are real and hence s n = \ + i^y n , for 7„ > T Q . 

Starting with a judiciously chosen Jost function, M(s, k), which is constructed 
such that M(s, 0) is Riemann's £(s) function, we have used inverse scattering meth- 
ods to actually construct a U(s; x) with the above properties. By necessity we had 
to generalize inverse methods to deal with complex potentials and a non-unitary 
S-matrix. This we have done at least for the special cases under consideration. 

For our specific example, f o °° \f {x)\ 2 V\(x)dx = 0, and hence we get no restric- 
tion on Img n or Res n . The reasons for the vanishing of the above integral are 
given, and they give us hints on what one needs to proceed further. The problem 
of dealing with small but non-zero energies is also discussed. 



1 



I. Introduction 

Many physicists have been intrigued by the Riemann conjecture on the zeros of 
the zeta function. The main reason for this is the realization that the validity of the 
hypothesis could be established if one finds a self-adjoint operator whose eigenvalues 
are the imaginary parts of the non-trivial zeros. The hope is that this operator 
could be the Hamiltonian for some quantum mechanical system. Results by Dyson 1 , 
and Montgomery 2 first made the situation more promising. The pair distribution 
between neighboring zeros seemed to agree with that obtained for the eigenvalues 
of a large random hermitian matrix. But later numerical work showed correlations 
between distant spacings do not agree with those of a random Hermitian matrix. The 
search for such a Hamiltonian in physical problems has eluded all efforts. Berry 3 has 
suggested the desired Hamiltonian could result from quantizing some chaotic system 
without time reversal symmetry. This seems to be in better agreement with numerical 
work on the correlations of the Riemann zeros, but one is still far from even a model 
or example. It is useful to explore new ideas. 

Our choice for this paper is an idea, due to Chadan 4 . In this approach one tries to 
relate the zeros of the Riemann zeta function to the " coupling constant spectrum" of 
the zero energy, ,S-wave, scattering problem for repulsive potentials. We sketch this 
idea briefly. 
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The Schrodinger equation on x G [0, oo) is 



d 2 f 



(A; k; x) + XV{x)f(X; k; x) = k 2 f(\; k; x), 



(1.1) 



dx 2 



where k is the wave number, A a parameter physicists call the coupling constant, V(x) 
is a real potential satisfying an integrability condition as in Eq. (2.2) below, and / 
is the Jost solution determined by a boundary condition at infinity, (e~ lkx f) — > 1 as 
x — > +oo. The Jost function, M(A; fc), is defined by lim x ^ f(X; k;x) = M(A; fc). It 
is well known that M is also the Fredholm determinant of the Lippmann-Schwinger 
scattering integral equation for S-waves. Both /(A; fc; x) and M(A; k) are for any fixed 
x > 0, analytic in the product of the half plane, Irak > 0, and any large bounded 
region in the A plane. In fact it is known that for any fixed k, Imk > 0, M(A; k) is 
entire in A and of finite order. Thus M(A; k) has an infinite number of zeros, A n (/c), 
with X n (k) — > oo as n — > oo. 

Starting with Eq. (1.1), and its complex conjugate with k = ir, r > 0, and setting 
A = A n (ir) we obtain 



For the class of potentials we deal with V = 0(e rnx ) as x — > oo. Thus we can take 



[7mA n (ir)] / |/(A„(ir); ir; x)\ 2 V(x)dx = 



(1.2) 



Jo 
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the limit r — > 0, and get 

/•oo 

[lm\ n (0)\ \f(X n (0);0;x)\ 2 V(x)dx = 0. (1.3) 

Jo 

Hence for repulsive potentials, V(x) > 0, all the zeros A n (0) are real. For any r, r > 0, 
the same is true for all \ n (ir). But \ n (ir) must be negative since the potential 
[A n (iT)V] will have a bound state at E — — r 2 , and that could not happen if V > 
and A n (ir) > 0. Hence by continuity, A n (0), for all n, is real and negative. 5 The zero 
energy coupling constant spectrum, A„(0), lies on the negative real line for V > 0. 
Chadan's idea is very simple. He introduces a new variable^, and defines 

A = s(s-1). (1.4) 

Thus one can write 

M(A,0)=M(s(s-l);0) = x(s). (1-5) 
It is easy to see now that for \Ims\ > 1, the zeroes, s n , of x( s ) are & h suc h that 

Sn = 7, + hn] A n (0) = S n (s n - 1). (1.6) 

The problem is actually somewhat simplified by noting that first we do not need 
the condition A n < as long as we restrict ourselves to the strip < Res < 1, and 
Ims > 1. Second, it is sufficient to prove that the integral in Eq. (1.3) does not 
vanish. Thus, one does not need a fully repulsive potential for the Riemann problem. 
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One might comment that it is very difficult to find a potential with A n (0) = 
s n (s n — 1) and s n — \ ± i~f n , s n being the Riemann zeros. But it is probably as 
difficult as finding a hermitian operator whose eigenvalues are j n . Indeed the latter 
may be impossible without introducing chaotic systems. 

The results mentioned above also apply when V — V Q + AVi, with only V\ > 0, 
and Vo, V\ both real and satisfying Eq. (2.2) and with certain restrictions on V Q . This 
remark leads directly to the basic idea of this paper. 

The objective of this paper is to show that the coupling constant approach can 
be significantly simplified and made amenable to inverse scattering methods. 

Our first remark is that one does not need a potential, V = V Q + AVi, depending 
linearly on the coupling parameter A. Given a one parameter family of complex 
potentials, U(s;x),x G [0, oo), which for fixed x are analytic in s in the strip, < 
Res < 1, Ims > T Q > 2, we can following similar arguments as above, obtain, for 
s = s n , s n being a zero of the zero energy Jost function, 

J \f(s n ,0;x)\ 2 [ImU(s n ,x)]dx = 0, (1.7) 

where / is the zero energy Jost solution evaluated at s = s n . 

Next, suppose in addition to the above properties, U has an asymptotic expansion 
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in inverse powers of s, actually better, s(s — 1), i.e. 



U(s, x) = V (x) + gVxix) + g 2 V 2 (x) + ... + g N V^(g; x); 



(1.8) 



where 



1 



(1.9) 



9 = 



8(s -iy 



Under suitable conditions on the V n (x) and estimates of the 0(g N ) term, and its 
phase, one gets again, 



with g n = [s n (s n — l)] -1 , the s n 's are the zeros of M(s, 0) the zero energy Jost 
function, and f(g; k; x) is the Jost solution with the full U . The result (1.10) is only 
established for zeros with Ims n > T Q , where T Q is large enough for the V± contribution 
to Eq. (1.7) to dominate the integral in Eq. (1.7). However, this is sufficient since 
the Riemann hypothesis has already been proved for zeros with \Ims n \ < T, where 
T could be as large as 10 5 . 

Again, all we need for, s n = \ + i^ ni is to have the integral in (1.10) not vanishing. 
In the end only the properties of V\ matter. 

In this paper we will use inverse scattering methods, albeit for complex potentials, 
to actually prove the existence of such a U(s; x). This potential, by construction, has 




(1.10) 
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the additional property that the zero energy Jost function is Riemann's £ function, 

\imM(s;k) = 2C(s). (1.11) 

fe^O 

We will also give explicit expressions for V , Vi, V 2 , and bounds on 
The difficult point turns out to be that in our specific example, 

roc 

/ |/(0;0;a;)| 2 Vi(a;)cfa; = 0. (1.12) 

Jo 

Thus we get no information on [Img n ], or Re(s n — \). We shall discuss what one 
needs to proceed further. This will require working with small, but non-zero energy 
values. 

We start by introducing a special class Jost functions, M ± , which depend on an ex- 
tra parameter v = s— \, with the property that the zero energy limit, limk^oM ± (u, k) = 
2£(f + |). For fixed v, the Jost functions are taken to be of the Martin 8 type, i.e. 
having cut plane analyticity in the momentum variable k. This is the class of Jost 
functions that results when the potential is a superposition of Yukawa potentials. We 
then use inverse scattering methods to prove the existence of a complex potential 
U (u, r) which is determined uniquely by the initial S'-matrix. We do carry out the 
analysis for v in the truncated critical strip, i.e. — \ < Rev < |, and Imv > T Q , with 
T > 16tt 2 . This is of course the domain most relevant to the Riemann problem. Stan- 
dard techniques of inverse scattering are not immediately applicable, because S(u, k) 
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does not satisfy the reality condition, and is not unitary for complex v. However, 
we shall see that in our specific case we can bypass these difficulties and carry out 
an inverse scattering procedure anyway. We have attempted to make the paper self 
contained, and do not rely on results that need the unitarity of S in the proof. 

In Section II, we give a brief review of relevant scattering theory results intended 
for mathematicians not familiar with them. This review also helps define our physics 
terminology. 

Section III is devoted to the introduction of our special class of Jost functions, 
M ± (u;k). Following that, in Section IV, we briefly discuss the real v case, which 
is a standard inverse scattering case covered by well known results. This section is 
instructive, even though real v is uninteresting for the Riemann problem. The next 
step, section V, is to study in more detail the properties of M ± . The main result is 
an asymptotic expansion in powers of a variable, g = (u 2 — which gives 

M (±) = M^(k) + gM[ ±) (k) + g 2 M^ + - + g^B^ig; k) (1.13) 

Here all the can be computed exactly via recursion formulae, and in addition, 

they satisfy [M^(k)\* = Mj~\k), and M^(-k) = Mj~\k). The remainder func- 
tions, are given explicitly and are 0(g) as g — > 0. 

In Section VI, with fixed v in the strip, we determine the number and positions 
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of zeros in the upper half fc-plane. It turns out that there is at most one such zero 
and it lies close to the origin. In fact we can give a good estimate of its position. 

Section VII is devoted to the study of the case \u\ — > oo, i.e. \g\ — > 0. Here 
M^(u,k) -> MW(fc) which is a known rational function in k. This leads to an 
exactly soluble Marchenko equation and an exact result for the corresponding V (x). 

Section VIII is devoted to proving the existence of solutions of the Marchenko 
Equation for our specific class of S'-matrices. With the resulting Marchenko operator, 
A(u; x, y), which is now complex, we proceed to define in the standard way a potential 
U(v,x) and corresponding Jost solutions, f^(u;k,x), of the Schrodinger equation. 
Finally, we check directly that indeed are solutions of the Schrodinger equation 
with the desired asymptotic properties. The main difference from the standard case 
is that U (v; x) is now complex unless v is purely imaginary. 

In Section IX we discuss the case v — it, t real. This is a standard inverse problem 
with S(it, k) unitary for k real, and the resulting U (it; x) is real. 

More detailed properties of V(v,x) are given in Section X. There we give an 
asymptotic expansion, V(v,x) = V (x) + gVi(x) + g 2 V 2 (x) + ... + g N vj^\g,x) with 
all V^s real and all representable by superpositions of Yukawa potentials. Also V n (x) 
is continuous and differentiable for xe[0, oo), and V n (0) is finite. For completeness we 
calculate Vi(x) explicitly, and indicate how V n (x), n > 1, can easily be computed. 
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We also give some needed properties of Vr and of [ImV^ ] for small Rev. 

In Section XI, we study the zeros, v n (k), of M~(u, k) for small fixed k with Irak > 
0. We prove that u n (0) are the standard Riemann zeros, and also that \u n (k)— f n (0)| = 
0(k 1 ^ p ) for small k. Here p is the multiplicity of the Riemann zero v n = u n (0). We 
also prove that any Riemann zero, Vj, is the limit of a zero of M^~\v\ k), Vj(k), as 
k -> 0. 

Finally, in Section XII we discuss the relation of our potential, V(g; x), and its Jost 
solutions to the Riemann hypothesis. We prove that in this case / °° |/(0; 0; x)\ 2 V\(x)dx = 
0, and hence no information on the Riemann hypothesis can result directly from this 
example at zero energy. But the reasons for the failure are clear, and they indicate 
the properties of a desired Jost function that will be sufficient to make the important 
step. The fact that one can set k = it, r > but small, and try to prove the hypoth- 
esis for u n (r), t arbitrarily small, but r ^ 0, provides a significant simplification of 
the problem. 
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II. A Sketch of Scattering Theory 



This section is intended to facilitate the reading of this paper by those mathe- 
maticians (or physicists) who are not familiar with elementary scattering theory in 
quantum mechanics. At the end of this section we will give a list of books and review 
papers where more information can be obtained. 

The Schrodinger equation for s-waves is given by 



+ gV(x)f = k 2 f, k = K + tr. (2.1) 



dx 2 

Here x G [0, oo), V(x) is real, g is a parameter that physicists call a coupling 
constant. The reason for introducing it will become apparent below. One studies the 
class of real potentials, V(x), which are locally summable functions and satisfy the 
condition, 

poo 

/ x\V{x)\e ax dx = C < oo, < a < m. (2.2) 

Jo 

For scattering theory the important solutions of Eq. (2.1) are the so called Jost 
solutions. 9 . These are the two linearly independent solutions, f^(g, k, x) with bound- 
ary values at infinity given by 

Jim e ±ikx f {±) = 1 (2.3) 
11 



Using the method of variation of parameters, we can replace Eq. (2.1) and the 
condition (2.3) by an integral equation: 

f^(g; k; x) = e* k * + g jT ™ k & - X) V(x')f^(g; k; x')dx> . (2.4) 

Starting with the papers of Jost 9 and Levinson 10 , the existence of solutions to Eq. 
(2.4) and their properties have been well established for V(x) satisfying the condition 
(2.2). 

The basic input needed is the upper bound on the kernel, 
sin k(x' — x) 



\r\\x'-x\ 

< Ci rrr, r , Irak = r, (2.5) 

- l+kx'-xV ' V ; 



k(x' — x) 

where C\ is 0(1). With this bound and the bound (2.2) one proves the absolute 
convergence of the iterative series of the Volterra equation (2.4) for any x > 0, and k 
with Irak > — y for f(~\ and Irak < y for f^ + \ Also it is easy to prove that for any 
finite g and x > 0, f( + \g; k;x) is an analytic function of k for Irak < y. Similarly, 
f(~\g; k;x) is analytic in Irak > — y. In addition, for k in the analyticity domain, 
the power series in g obtained by iterating Eq. (2.4) is absolutely and uniformly 
convergent for g inside any finite region in the g-plane. Thus both f ± (g;k;x) are 
entire functions of g. 

The scattering information is all contained in the Jost functions, denoted by 
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M^(g; k) and denned by, 

M^(g;k) = \imf^(g;k;x). (2.6) 

x— >0 

Both limits in Eq. (2.6) exist for finite \g\, and in the respective domain of analyt- 
icity, for all potentials satisfying the condition (2.2). The S-matrix is given by 

^'V- M(-)(g;kY {ZJ) 
For real g and Imk > 0, M^(g, k) has no zeros except for at most a finite number 
on the imaginary k-axis. These zeros, k n = ir n , give the point spectrum of the 
Hamiltonian of (2.1) with E n = — t%. Their number cannot exceed the value of the 
integral / °° x\ V\dx, a result due to Bargmann 11 . 

Another important property of M^~\g\ k) was first obtained by Jost and Pais 12 . 
The regular solution of Eq. (2.1), (f)(g; k; x), with (J)(g; k; 0) = 0, is 

<j>(g; k; x) = 7^[M {+ \g; k)f -\g; k; x) - M^(g; k)f^(g; k; x)}. (2.8) 

The solution <p satisfies a Fredholm type integral equation which, for potentials 
satisfying (2.2), was studied in ref.12. Jost and Pais demonstrate explicitly that 
M(~\g, k) is identical to the Fredholm determinant of the scattering integral equation 
for <fi. Hence for any fixed k, with Imk > 0, the zeros of M^~\g\ k) in the g-plane, 
g n {k), give the "coupling constant eigenvalues" at which the homogeneous Fredholm 
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equation has solutions, = g n {k)Kcf). Since, M^~\g, k) is an entire function of finite 
order in g, the sequence gi(k),g 2 (k), ...,g n (k) tends to infinity as n — > oo. 

For the purposes of this paper a result of Meetz 5 is instructive. Let us consider 
a potential which is repulsive, i.e. V > for all x G [0, oo). Then for k = ir,r > 
the coupling constant spectrum, g n {ir), is real and negative. This result is implicitly 
contained also in ref. 12. 

In this brief review we need to make an important remark about complex po- 
tentials, V ^ V* . Mathematicians and mathematical physicists often ignore these 
potentials. The Hamiltonian is no longer self-adjoint if V ^ V*, with g — 1. But 
physicists, especially those who work on nuclear physics, do not have such a luxury. 
There are many interesting and useful models, especially in nuclear physics, where V 
is complex. Of course, the general and beautiful results which hold for real V do not 
all apply for complex V. But many survive, and one has just to be careful which to 
use and to establish alternative ones when needed. 

There are many books that cover inverse scattering. But for the purposes of this 
paper we recommend the book of Chadan and Sabatier 13 , since it also discussed the 
superposition of Yukawa case and the Martin results. For the standard results on 
inverse scattering the review paper by Faddeev 14 is highly recommended. 
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III. A Special Class of Jost Functions 

In this section we will combine two results whose progeny could not be more 
different to obtain a representation for a class of Jost functions that we shall study 
in detail. The first is Martin's representation for the Jost functions of the class of 
potentials that can be represented as a Laplace transform. The second is Riemann's 
formula for the function £(s) defined below. 

Starting fourty years ago, physicists 15 ' 16 , for reasons not relevant to this paper, 
studied the class of potentials that in addition to satisfying Eq. (2.2) have a Laplace 
transform representation, i.e. for all x > 0, 

V(x) = / C(a)e~ ax da, m > 0, (3.1) 

J m 

where C(a) is summable and restricted to satisfy J~ \C(a)\a~ 2 da < oo. This last 
condition guarantees that x|V"(a;)| is integrable at x — 0. 

For these potentials Martin 17 proved that the Jost functions M^(k) have the 
representation 

M( ± )(Jb) = l+ r^Tda. (3.2) 
if a±ik 

Here w is real and summable and is such that M^> — > 1 as \k\ — > oo. We have 
set g = 1 here. Note that not any arbitrarily chosen summable w(a) is acceptable. 
M^~\k) must have no zeros for Imk > except for a finite number on the imaginary 
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/c-axis corresponding to the point spectrum. 

For our purposes here we choose a specific family of functions M^(u; k) defined 
such that 

M (±) (^0)e2((,^), (3.3) 

where 

m = \s{s - l)7r-f r(|)C(s) (3.4) 

and 

oo 

((s) = Y / n~ s ] Res>l. (3.5) 

n=l 

Riemann's formula for £(s) defines an entire function of order one in s, and is 
given by 

2f(s) = l + s(s- I) J™ i){a)[a^ 1 + a^'i}da, (3.6) 

where 

oo 

^(a) = e~ nn2a ; a > 1. (3.7) 

n=l 

We also have the symmetry relation £(s) = £(1 — s). 
For convenience we define the variable, v, as 

s = \ + v (3.8) 
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With this variable + v) is symmetric in z/, and we have 

2f(- + v) = 1 + [v 2 - -) / ^(a)a"J[a^ + a^]da. (3.9) 

Zi 4r t/ 1 

Our starting point is to define two functions, M ± (u; k), as 

M<±V; A;) = 1 + (z/ 2 - \) f jKgjglM+ajr] ^ (3 . 10) 

This definition holds for any finite, real or complex, is, and for any k excluding the 
cuts on the imaginary A;-axis, k = it, 1 < r < oo, for M^, and — 1 > r > — oo, for 
MM. 

Obviously, we have 

M( ± \v;0) = 2a\ + v). (3.11) 
In addition, the fact that ip(ct) = 0(e~ na ) as a — > +oo, guarantees that for any finite 

IH 

lim M ± (u;k) = 1. (3.12) 

This is true along any direction in the complex plane excluding the pure imaginary 
lines. But even for argk = ±| the limit holds using standard results. The immediate 
question that faces us at this stage is: for which regions in the //-plane, if any, can 
one use the functions M^(u; k) defined in Eq. (3.10) as Jost functions and proceed 
to use the resulting S-matrix, S(u; k), as the input in an inverse scattering program. 
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There are two issues involved. The first, and most important, is to make sure that 
M^~\v\ k), has no complex zeros in k for Imk > 0, except for a finite number on the 
imaginary axis. This is not true for any v. But fortunately for the set of z/'s most 
important to the Riemann hypothesis M^~\v\ k) has at most one zero close to the 
origin with Imk > 0. This will be shown in Section VI. 

The second issue relates to the question of reality. For real potentials V and real 
k we have the relations [M^ik)}* = M^(k), and \S(k) \ = I. Clearly for complex v 
this does not hold for M^\v\ k). However, we will prove that for those values of v 
in the truncated critical strip, one can still carry out the inverse scattering program 
and obtain a unique and well defined V{y\ x) which of course now could be complex. 
Since the old results of inverse scattering theory all use the fact that \S(k)\ = 1, we 
have to go back to square one and prove every step anew for the present case. Our 
task is tremendously simplified by the fact that even though S(v; k) is not unitary, 
we still have \S(u; k)\ = 1 + 0(l/\u\ 2 ), and we are only interested in \u\ > 10 3 . 
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IV. The Real v Case 

For potentials satisfying the representation (3.1), Martin 19 , in addition to the 
results summarized in Eq. (3.2), developed an iterative scheme which enables one to 
reconstruct the measure, C(a) in Eq.(3.1) from the knowledge of the discontinuity 
of S(k) along the branch cut on the imaginary axis, k = ir, y < r < oo. This 
gives an inverse scattering method that at first sight looks quite different from the 
standard ones of Gelfand-Levitan and Marchenko. The relation between these two 
methods was first clarified by Gross and Kayser 20 and independently by Cornille 21 . 
They showed that for potentials of the form (3.1) the Marchenko kernel is a Laplace 
transform of the discontinuity of S(k), and they carried out an extensive analysis of 
the relation between Martin's and Marchenko's method. These results were reviewed 
and enlarged in a more recent paper by the author 22 . 

For v real and \v\ > |, the functions M^(u; k) defined in Eq. (3.10) are indeed 
bona fide Jost functions with (M^ + \u; k))* = M^~\v\ k) for real k. The positivity of 
i]){a) guarantees the absence of a point spectrum. The S-matrix is 

We define the discontinuity, D(u; r) as 



D(u,t) = \im[S(u;iT + e) - S(u;iT - e)}, r > 1; (4.2) 
19 



From Eq. (3.10) one obtains 

D(u, r) = ! J{U ' ,T ) m , (4.3) 

with 

uj(v,t) = 7r(z/ 2 - 1/A)i/j(t)t^[t^ + T~]. (4.4) 

For real v > |, _D(z/, r) and u(v, t) are positive for all r > 1. The case with u;(i/, r) > 
is the easiest to handle by the Martin inverse method, and it can be done explicitly. 

Although having v real and v > | is of little direct interest to the Riemann 
problem, we give the results here as they might be helpful to the reader. For details 
one should consult ref. 22. 

The S-matrix, S(u,k), uniquely determines a potential, V(u;x), and its Jost so- 
lutions f^(iy; k;x). V is given by 



00 —1 r°° r°° (Y\ n n D(v a )e~ 2ajX ) n 

v(*)=*u-r»i *>..../ ^ "fe ( a; + i w )) ( £ aj) - (4 ' 5) 

This series for V is absolutely and uniformly convergent for all a; > 0, and v > |. 
This follows from the positivity in (4.3), 

1 /■«> \D(v,a)\ , ^f^°u(u,a)/ada , t , 

7l cfa < s J1 1 v /' < 1. (4.6) 



tt7i a + r - i + I « d/j 
The Jost solutions, j^{y\ are given by 

/(± , = ^ + eWte £ ( zi r . r dQ „... r da „ r J^^f:**?.,, (4.7) 



re =o n Jl Jl [U.]=o(aj + a j+1 )}[a ±ik] 
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Again this last series is absolutely and uniformly convergent for all x > 0, and A; in a 
compact domain inside the respective regions of analyticity. 

One can check directly that f^ 1 given by Eq. (4.6) are solutions of the Schrodinger 
equation with V{y\ x) of Eq. (4.4) as potential, see ref. 22 for more details. 
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V. Some Properties of M^\u, k) for |z/| > 1. 

To proceed further and study M^(v,k) for complex is, and more specifically v 
in the critical strip,— \ < Rev < \\ Imv > 1, the defining representation (3.10) is not 
fully instructive. This is because the behavior of M ± for large Imv is not adequately 
shown by Eq. (3.10). Our final result in this section is to obtain an asymptotic 
expansion of M ± (v, k) for fixed k in inverse powers of [y 2 — |). 

We need to carry out integrations by parts on the integrand in Eq. (3.10) analo- 
gous to those performed in Titchmarsh's book 18 , for Eq. (3.9). 

The following lemma will prove extremely useful. 

Lemma 5.1 

Let W(a),a G [l,oo), be a C°° function, and W(a) = 0(e~ a ) as a — > oo, then 
given the integral 

W{a)[a% +a^}da, (5.1) 
one has after two integrations by parts 

I(v) = 2 _ {8[W(1) + (W'(a)) a=1 \ + / ^(a)[a5 + a^]da}, (5.2) 

where 

Wi(a) = ^W(a) + 12aW'(a) + 4a 2 W"(a). (5.3) 

Proof: 
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We rewrite Eq. (5.1) as 

/OO 3 v 3 — v 3 

da(VF(ai)a!4)[a!2 - 4 + a — - 3 j 

Integrating by parts we get, 

2W(1) [°° ( d 3 lu ai + i a^+i 

*M - /„ 2 _ i\ ~ / ^(T-^^l^DlT-r-Trr]' ( 5 - 5 ) 

V 4/ tt " 2 ~ l ~ 4 2 4 

This again can be rewritten as 



-v 5 



w . 2W(1) f 00 , , 3,6? rTTr , . 3 1N1r Q;2 4 Q;— 4 

^ = T^rtk ~ \ da { a H^-[W(a)a-*])}[ 1 — T - 17 — H. (5.6) 
^ 4 ; ji aa 2 4 2 4 

Carrying out a second integration by parts we obtain, 

= ^r[2^(l) + 8{^(Vy(a)af)} a=1 ] (5.7) 

Performing the differentiations in (5.7) easily leads to Eq. (5.2). 

We can apply this lemma to the integral in Eq. (3.10) which defines M^\v, k). 
Setting 

W^^, (,S) 
and restricting k to the corresponding domain of analyticity in k 

P (+) = {k\Imk < 1}; P ( - } = {k\Imk > -1}. (5.9) 
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we get 

M (±) (z/; k) = M^\k) + daW[ ±] (a- k)[a% + n ~\. 



with 



1 " {^(a±ik) v 



and 



The first term in (5.10) is independent of v, and given by 

ai a 2 



M<±>(*) = ! + -£- + 



l±ik (l±ik) 2 ' 
with 

ai = -1+8^(1) ; a 2 = -8^(1). 

In obtaining (5.14) we have used the identity 

4^(1) +^(1) = -1/2. 

It is important to note that both a± and a 2 are negative and that a x + a 2 
leads to 

limM (±) (£;) = 0. 
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As a check on Eq.(5.10) we take the k — > limit, 

#'(i/;0)= rdaE^/a^laHa^] (5.17) 

Substituting the expressions for &p given in Eq. (5.12) we get 

M (±) (z/,0) =4 / da(il>"(a)a* + -ip'(a)a*)[a% + a^} =2£(- + v). (5.18) 
Ji 2 2 

The lemma 5.1 can be used repeatedly to give an asymptotic expansion of M ± (u; k) 
in inverse powers of (z/ 2 — |). Recursion formulae can be given to give each term from 
the preceding one. 

Indeed given 

fc) = p - i)n-i / ^ ±} («; *)[«' + «^], (5.19) 

with 



2n+l J,!")/ ^ 



and G P ± , one obtains 



45 = . 2 1 1)W {8[WW(1; fc ) + (^^) a =i]+/°° dQWSUa; fc)[aUa^]}, (5.21) 



(^ 2 - I) 

and 



W£l = ^W™ + ^ + 4a 2 (^))"; (5.22) 
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where the primes indicate differentiation with respect to a. Again W^+i) w ^ be as 
in Eq. (5.20) 

2n+3 $(; t + 1 )f a ) 

K + » = g (5-23) 

For each n we have 1 < £ < 2n + 1, and the functions $^ satisfy a recursion 
formula, which follows from (5.22). 

$5 n+1) (a) = ^5 n) + 12a(^ n) )' - 12a(* - (5.24) 
+ 4a 2 [($5 n) )'' - 2(€ - 1)($S)' + (€-/)('- 2)<&&] 

All the can thus be determined by iteration starting from 

$S 0) (a) = ^(a)ai (5.25) 
The general form of ^\a) is easily determined to be 

2n+l-£ 

$W( a )= £ C (n) (^;j) alpha^+i-^ia). (5.26) 
i=o 

The coefficients C^ n \l;j) are real, and C^(l; 0) = 1 determines all the others. Also 

^\a) = (-^) j il>(a). (5.27) 
da 

At this point we can substitute Eq. (5.26) in (5.24) and obtain a recursion formula 
for C^(£;j), 
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C W(t,j) = C^(£;j)[^ + 12(±+£ + j-l) (5.28) 

- C<»>(€ - 1, j)[12(€ -l) + 8(€-l)(|+€ + j- 2)] 

+ C (n) (£, j - 1)[12 + 8(i + £ + j - 2)] + 4C< n >(4 j - 2) 

- 8(£-l)C {n) (£-l;j - l)+4(£-l)(£-2)C {n) (£-2;j). 

Here we have 

l<£<2n + l, <j <2n + l- £. (5.29) 

For all other values of I and j, C^(£,j) = 0. 
Starting with 

C<°)(1,0) = 1, (5.30) 

we can compute all other C( n )(£, j). For example,C«(l, 0) = 6, C«(l, 1) = 14, and 
(7^(1, 2) = 4. This agrees with the direct calculation given in Eq. (5.12). In table I, 
we give all the coefficients C^ n \£, j) up to n = 4. All the coefficients are integers. 
Finally, we give the general form of the surface term in Eq. (5.21). We define 

M^\k) = 8[^)(U) + (^)a=i]. (5.31) 
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From Eq. (5.20) and (5.26) we obtain after some algebra, 



2n+2 (n) 

= g p^, (5.32) 



with 



2n+2-^ I 

xi n) = 8 ]T CW(£,i){(- + £ + j)V' (j )(l)+V' (j ' +1) (l)} (5-33) 

3=0 4 
2n+2-^ 

- 8 J2 (t-i)& n \e-i,j)i><f\i). 

3=0 

For the purposes of this paper it is sufficient to apply our lemma up to the n = 3 
level. We introduce as a new variable, g, 

9 - ^ (5.34) 

Our final result for M (± (z/, k) with A;eP (±) and >> 1, 

AfW(i/; fc) = M^\k) + tfM^ife) + (? 2 M 2 (±) (A;) + g 2 B^\v, k). (5.35) 

Here we have 

4 6 



with be = and 



M 1 (± 'W=g[ r ^ p , (5.36) 



Mf»W=g lr ^ p , (5.37) 
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C£ = Xe ■ The remainder function R 2 is given by 

roc Z ( rv\ 

Rf\v-k)=J i (E|^p)[« f + <*?]da. (5.38) 

For \Imv\ > 10 3 , the first two terms of Eq. (5.35) give a very good estimate for 
M^l We shall explore this in much more detail later. One can go to higher orders 
in g, but the resulting series is only asymptotic. For our purposes here Eq. (5.35) is 
enough. 

It is important to stress another property of and namely as k — > 0, 

M 1 (±) (0) = 0, (5.39) 
M 2 (±) (0) = 0. 

We have already shown that M ^(0) = 0. To check this we give the explicit form of 
the coefficients be in (5.36). Using — be, Eq. (5.33), and table I, we get 

6i = 60^(1) + 300^(1) + 216^"(1) + 32^ w (l), (5.40) 

b 2 = -300^(1) -432^(1) -96^"(1), 

6 3 = 432^(1) + 192^(1), 

64 = -192^(1) 

Numerically, the fe's are given in table II. 
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We now have 



Ml (±) (0) = X> (5.41) 
1=1 



1 ^ 1 ^ 

32[^'(1) + ^"(1) + -^'(1)]. 



But 



+ ^V(l) + = (5.42) 

This identity follows from the relation 18 

y/a(2il)(a) + 1) = 2^(l/a) + 1. (5.43) 

Differentiating (5.40) once and setting a — 1 immediately gives Eq. (5.15). Differen- 
tiating three times leads to Eq. (5.42). 

Indeed there is an infinite sequence of identities like Eq.(5.42), always starting 
with ip 2n+1 (l), odd derivatives, which result from differentiating Eq. (5.40) (2n + 1) 
times. Thus, again 

M 2 (±) (0) = ^> = 0, (5.44) 

depends on the next identity: 

^)(1) + |^)(1) + ^V)(l) + ^W(l) + ^§V>(1) = 0. (5.45) 

Of course only the first two coefficients in (5.45) are unique, since we can always add 
a multiple of the l.h.s. of (5.42) to (5.45). 
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The vanishing of M^(0), j = 0, 1, 2, is indeed necessary since 

/i^V; 0) = 9 2 r da(t ^f 1 )^ + = n\ + »). (5-46) 

which is the result of carrying out four more differentiations by parts on the formula 
for £{\ + v) given on page 254 of reference 18. 

In table II, we give the numerical values of Xt™' \ f° r n — 1> 2, and 3, and b e = x^P 
while ce = x~P ■ 

For the convenience we summarize the results of this section: 

N 

M&(v, k) = £ rfW + ^^V, (5-47) 



ra=0 



where 



= (z/ 2 - 1/4)" 1 , (5.48) 



2n+2 (n) 



and are real numbers given in Eq. (5.33). In addition we have 

2n+2 

E X? = 0, (5.50) 
e=i 

which guarantees that M^(k) — > as — > 0. For real fc, we have [M^ + ^(A;)]* = 
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Finally the remainder term Flffi is given explicitly by, 



m v ,k> =/ ^g^J^K^^i, (5 . 51) 

with 

2n+l-£ 

$f } (a)= £ C (n) (£;j> 1/4+ ^-V (i) (a), (5-52) 

3=0 

The C^ n \i,j) are integers determined by a recursion formula given in Eq. (5.28), 
with C(°)(l, 0) = 1, and tp (j) (a) are the jth derivatives of i/j(a), Eq. (5.27). 
For k = 0, we have 

M (±) (z/, 0) = 2^(z/+ 1/2). (5.53) 
From Eqs. (5.47) and (5.51) we then have for any integer n > 0, 

1 2n+3 \ 

2£(^ +^)=9 n da( £ \^fV /2 + a""/ 2 ]. (5.54) 

For n = 0, this formula is given in ref.(18) on page 225. The results for larger n can 
be obtained by successive integrations by parts. 
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VI. The Zeroes of >(y \ k) for Imk > 0, and fixed v 



To study the Riemann hypothesis we need only to focus on the truncated critical 



Since the Riemann Hypothesis has already been rigorously established up to Imv = 
0(1O 6 ), we can simplify the calculations of this paper tremendously by taking T Q to 
be large. Initially, we take T Q = 10 3 . 

The following lemma will be quite useful. 

Lemma 6.1 

For any v G S(T a ), and k such that Imk > —\, we have 



strip, S(T ), 



S(T ) = {v\Imv > T c 




(6.1) 



\M(-\v,k)-M{,-\k)\ < 



C 2 



(6.2) 



and 



C 2 < 10 3 . 



(6.3) 



Proof: 



Taking the expansion of \v\ k) in powers of g — {y 2 



j) 1 to first order we 



have 



M ( -\u, k) - M^~\k) = gM[~\k) + gRt\u, k), 



(6.4) 
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where M{ (k) is given by Eq. (5.36), and 

5 (2) 

J 1 £ \CX % /CJ 

(2) 

Here is given by Eq. (5.26) and table I. 
First, we have for Imk > — |, 

\Mt\k)\ < £ JTrLl ^ E N-d)']- (6-6) 

Using table II, we get 

IM^A;)! < 68, Imk > ~. (6.7) 
The upper bound on R[ ^ for Imk > — \ is 

\R[-\u,k)\<2 aHE\*?\<*)\<*- e (a)')> ( 6 - 8 ) 
- 71 e=i 6 

where we have used \a/(a — ik)\ < | for a > 1 and Imk > — |. Using Eq. (5.26) we 
have 

ME(o)'l E \C^(£,j)\a^^\a)\} (6.9) 
where we note that C^ n \£,j) = (— l) f+1 .|C^(£, j)| as can be seen from Table I. 

The series in Eq. (3.7) that defines ip(a) is highly convergent for a > 1. Indeed 
the first term gives a good approximation to it and to its first six derivatives. One 
can easily derive the bounds, < j < 7, 

n j e- na < \i) {j \a)\ < 7r j e- na (l + e(j)), (6.10) 
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where e(j) is 



e(j) = e-*<[l + (2n 



(6.11) 



For j < 4, e(j) < 0.021. Thus it is sufficient for the purposes of this estimate to use 
ip(ri(a) = (—iy iri e~ na . Substituting this in (6.9) and carrying out the a integration 
we get 



where T(j, (3) is the incomplete gamma function. From Table I, it is now easy to check 
our bound of C 2 = C 2 + 68 < 200. This completes the proof of lemma 6.1. 

It should be apparent to the reader that one could use more refined methods to 
obtain a much better bound on r[~\ We do not do this at this stage. Our most 
important task is to study the Riemann conjecture for Imv > T Q with T Q taken below 
the maximum for which the hypothesis has been rigorously established. In a future 
paper, we will try to find the lowest value of T Q for which our method works. 

The function M^~\k) given by Eq. (5.13) is a rational function of k 



R[-\v,k)\< 




t=\ u j=0 



(6.12) 



Mt\k) 



1 + 



ai (-1 - ai) 



(6.13) 



l-ik (1-ik) 2 ' 
k^ + ii^ + a^} 



(l-ik) 2 



Here 



ai = 



1 + 8^(1) = -0.6543 



(6.14) 
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Obviously, has two zeros, k\ = 0, and k 2 = —i(2 + ai) = + 8-0(1)). Thus 

Imk 2 < —1. Hence, M^~\k) has only one zero in the half-plane, Irak > 0. 
Focusing on the domain Imk > —\, and \k\ > \ we get, with k = k + it, 



, . \Jk 2 + t 2 Jk 2 + (r + ri) 2 

i^-'wi = — K , + v + iV ■ (6 ' 15) 

where rj — 2 + a± > 1. 

It is easy to find a lower bound for \M^ | in the above domain. Setting q = (0, —77), 
and p = (0, —1) we get 

\Mt\k)\ = E^A > J^l_ > _L (6 .i 6) 
\k — p\ 2 \k — p\ 2 25 

We can now prove the following. 
Lemma 6.2 

For any v £ S(T Q ), and k such that Imk > 0, with \k\ > |, M^~\u,k) has no 
zeros and has a lower bound 

|M(~V,fc)| > (0.04-1). (6.17) 



- 



Proof: 



From lemma (6.1) we have 



|M<->(i/,fc)|>|M(->(*)|-^, (6.18) 



7^2 ' 

o 
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with C 2 /T < 2. Using (6.16) we get 

|M(-)(^)|>(0.04-1)>1 (6.19) 

Thus any zeros, k Q of M^~\u, k) in the upper half k-plane must have \k a \ < \. 
Proceeding further we have: 
Lemma 6.3: 

For any fixed v G S(T ), the maximum number of zeros ofM { ~\v] k)with Imk > 
is one. 
Proof: 

From lemma (6.1) we get for v e S(T ) 

|MH(^)-M(-)(^)|<1; \k\ = \, (6.20) 

By very similar arguments we can also show that 

.dM^(u,k) dM<r\ A ,,, 1 , 

where A = 0(1). 

Let ivj ^ denote the number of zeros of M^~\k) in the disc \k\ < j, and Ni be 

4 ' 4 

the corresponding number for M^~\u, k), ueS{T ), then 

N , - N r = ± i - M~)?h (6.22) 
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where Ci is the circle \k\ — 4, and the prime indicates differentiation with respect to 
k. 

Hence, 

|7V r M°)| < ^^{Max^^ i |M^ - ^ (A;) [M^ (z/, A;)]'— [M^ - - 1 (/c^'M^ (z/, A;) | } < (^) < 1, 

(6.23) 

where we have used Eqs. (6.16), (6.19), and (6.21). Since M^~\k) has only one zero 
in the disc so does M^~\v, k) and lemma (6.2) which proves the absence of zeros with 
Irak > 0, and \k\ > |, completes our proof. 

The next question is where is this one zero of M^~\v\ k)l To answer this question 
we first order expansion for M(-)(z/, k), 

s $f(a) 



M(")(i/, fc) = Mt\k) + ^"(fc) + g I " da(E + «^ ] (6-24) 



With /c = 0, we have 



2£(«/ + 1/2) = g r da(Z -*"P).[a* + cT?\ (6.25) 

*/ 1 £ ^ Qt 

Subtracting these two equations we get 

/'OO ^ 1 1 

M^\u,k) = 2i{u+l/2)+Ml\k)+gM[-\k)+g / da{Y,i , _ .u V --}*?\<*)[<**+<*^ 

(6.26) 

Note that now the integral on the r.h.s. is also 0{k) as k — > 0. 
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The £ function for large values of \Imv\ is exponentially small, mainly due to the 
T(|) factor in Eq. (3.4). In fact given standard results on the order of ((s) in the 
critical strip we have £(| + v) — 0(\u\ p e 'T") where p = 2 + 5 and < 5 < 1. 
Thus from Eq. (6.26) and the exact expression for M^~\k) in Eq. (6.13), we get the 
position of the zero, k , near the origin, 



where (2 + a x ) > I. Hence 



-2ite£(± + v) , 
/mfea (2 4- a.,) ' (6 ' 28 > 



If Re£,(\ + v) > 0, M^~\v,k) has no zeros for Irak > 0. On the other hand if 
v G S(T ) is such that Re^{\ + v) < there will be one zero close to the origin, but 
in the upper half k plane. 

Finally, we note one important fact, namely, if (u + l/2) is a zero of the {-function, 
then M(~\v , k) has no zeroes for Irak > 0, and its only zero with Irak > —e occurs 
exactly at k — 0. 

In summary M^~\v;k) , with ueS(T ), has most of the properties of a Martin 
type Jost function with the exception of one, i.e. reality. We list these properties: 

i) M^~\v, k) analytic in the cut k-plane with a cut for k = it, — oo < r < — 1. 

ii) \im lkHoo M(~\v,k) = 1. 
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iii) M( \is, k) has no zeros for Imk > |, and veS(%). 

iv) If we write 



I POO , _ 

u) = [v 2 - -) / dail>(a)a* [a* + a~ }e~ au , (6.29) 
4 ii 



then 



roo 

M {±) (v,k) = l+ Z(v,u)e Ttku du (6.30) 

Jo 

In fact we could have used Eqs. (6.29) and (6.30) as the starting definitions of 
Af( ± )(i/, fc). 

v) We define the S-matrix 
and it follows that 

/+oo 
\S{u,k) - l\ 2 dk < oo. (6.32) 
-oo 

vi) The reality condition does not hold for all v. If v is purely imaginary, i.e. v = it, 
then for real k, (M (+) (z/, k))* = M (_) (z/, fc), and hence \S(v, k)\ — 1. However, if z/ is 
non-, i.e. i?ez/ 7^ 0, then the above relation does not hold. However, we still have for 
real k 

\S(u,k)\ = l + 0(±). (6.33) 
For an arbitrary v, ueS(T ), we can still carry out the inverse scattering program 
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by properly handling the one zero in the upper half k plane. We will do that in 
Section VIII, where the resulting potential is complex. 



41 



VII. The Limit Case, |z/| — >• oo 

Before we proceed to the main proof, we shall solve exactly the limiting case 
\u\ — > oo. This result will be extremely useful in the rest of this paper. 
We start with 

M^\u,k)M^\k), |i/|->oo. (7.1) 
where from Eq. (6.13) we have 

M« W = =S^|±M . (7. 2) 

We thus have Jost functions which are rational in k. This is the case first studied 
by Bargmann 23 in the paper which gave the famous phase equivalent potentials. 
The S-matrix is also rational, 

S m = fl_ k-i{2 + ai ) l-ik 

where (2 + ai) > 1, ai = —1 + 8^(1). One can use Bargmann's method to uniquely 
determine a potential V (x) which has the S-matrix given here. But we prefer to 
determine V a by using Marchenko's method. 
The Marchenko kernel F Q is 

Fo(x) = — / (S (k) - l)e ikx dk. (7.4) 

This Fourier transform converges in the mean, (S — 1) — > 0(1/ k) as — > ±oo. 
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By contour integration 

F (x) = \ e~ x + X 1 xe~ x , (7.5) 

where 



^ = 8a ' + + 4 "\7 4 . 

(3 + ai) 2 



and 



(3 + ai) 

Both A D and Ai are negative. 
The Marchenko equation is, 

A (x,y)=F (x + y)+ A (x,u)F (u + y)du; (7.8) 

and A (x,y) = for y < x. With F G as defined by Eq. (7.5), one can easily obtain 
the exact solution of the integral equation (7.8). 

From Bargmann's paper 23 it is clear that we have the ansatz, 

A (x,y) = \B{x) + (y - x)C (x)]e~<*-*) . (7.9) 

Substituting this trial solution in Eq. (7.8) and carrying out the integration over u 
we get 

A (x,y) = e -0'-*){ e -^[(^M + ^l)(A o + A 1 x) (7.10) 

(B(x) + C{x)) ^-t. Mx)^ C(x)X 1 
+ ^ Ai + A + xAiJ + ye [ — 1 h AiJ}. 
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But from Eq. (7.9) we also have 



A (x,y) = [(B-xC)+yC]e-( y - x \ (7.11) 
Comparing the coefficients of y in Eqs. (7.10) and (7.11) we get 

e-CM-^ + ^i + A,, (7.12) 



and the terms to zero order in y give 

(f + j)(A + X lX ) + ^ + ^ + X + xX 1 = e 2x [B - xC\. (7.13) 
These last two equations determine B(x) and C(x) giving, 

C {x ) = AJ A '-< A ^'" , (7.14) 

and 



B(x) = A ; +2A ^ +( ^ , (7.15) 

Since X ± < 0, X Q < 0, and (Af/16) << 1, the denominators in Eqs. (7.14) and (7.15) 
do not vanish for any x > 0. 

One can simplify Eqs. (7.14) and (7.15) by defining 

p=-± = (l±^) = 0.1474. (7.16) 

H 4 K 3 + a 1 J V ; 



Then its is easy to show that 

1 



2 (A + A 1 ) = l-p 2 . (7.17) 
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We obtain 



C(x) = -(4p + Ve- 2 -) 

{ ) [e 2 - - p 2 e~ 2 - + (1 - p 2 ) + 4px] ' { } 

and 

s( , = (2p 2 + 4p-2)-8px + 4p 2 e- 2 - 
W [e 2 - - p 2 e- 2 - + (1 - p 2 ) + 4px] • 1 ' 

From Eq. (7.9) we see that the potential, V (x), is 

Vc{x)= _ 2 ^ (x , x)= - 2 tm. (7 , 0) 

This leads to 

(4 - 2p 2 - 16p)e 2x - 4(1 - p 2 )p 2 e- 2x - 16px(e 2x + p 2 e~ 2x ) - 64p 2 x - 16p 2 



V {x) 



[e 2x - p 2 e' 2x + (1 - p 2 ) + 4paf 



(7.21) 



Note that V Q (x) = 0(e~ 2x ) as x -> +oo. 
The two Jost solutions are given by 

fo (k,x) = e^ kx + / dyA,(x,y)e^. (7.22) 

Substituting Eq. (7.9) for A c , we get 

One can now check directly that 

+ V (x)ff = k 2 f±, (7.24) 
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also /W -> Mg(k) asx^O. 

The results of this section can also be obtained by using the technique developed 
by Bargmann 23 which preceded the results of refs. 6 and 7. The Jost functions 
defined in Eq. (7.2) do indeed determine uniquely the potential V (x) and its solutions 
/ (±) (£;, x). The fact that M^\k) = 0(k) as k -> 0, leads to the degenerate case using 
Bargmann's method but the final results agree with those by Marchenko's method 24 . 

Finally, it should be remarked that the full scattering amplitude for V (x) has a 
pole at k — 0. However, this is not part of the point spectrum. There is no L 2 (0, oo) 
solution of the Schrodinger equation with V (x) for k — 0. 
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VIII. The Marchenko Equation 

We define an S'-matrix as in Eq. (6.34) 

_ ffM 

MH(^)' (8 - lj 

where z/ G S(T ). S(v,k) is, for any u, analytic in k in the strip — 1 < 7mA; < +1. 
Next we define a Marchenko kernel for S(v, k), 

F(is; x) = f dk[S(u, k) - l]e lkx , x > 0, (8.2) 

where L is a line Imk = <5>0,1>5>0, and without loss of generality we fix 
5,5 = 1/4. This Fourier transform is convergent in the mean, (S — 1) = 0(1/ k) as 
Rek — > ±oo. 

Actually we can also perform an integration by parts on (8.2) for any x > e > 0. 
Since dS/dk is bounded and (dS/dk) = 0(1 /k 2 ) as k — > ±oo, we have absolute 
convergence for any x > 0. 

It is important to note here that as shown in Eq. (6.32), [S(u, k) — l] G L 2 (— oo, +oo), 
along the line Imk = 1/4. 

Of course Eq. (8.2) is not the standard definition of the Marchenko kernel. In the 
standard case one integrates along the real fc-axis. If we move the contour in (8.2) to 
the real axis, then there could be an extra contribution from the pole produced by 
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the zero of M~\u, k) when Re^(u + 1/2) < 0. But for a Martin type S-matrix all the 
scattering data, including that coming from the point spectrum, is contained in the 
discontinuity across the branch cut on the imaginary k-axis, see references 20,21, and 
22. 

We prove the following lemma: 
Lemma 8.1: 

F(v,x) is a.) continuous and differentiable in x, x G [0, oo); b.) F(v,x) = 
0(e~ x )as x — > +oo; c.) Both F(v, 0) and F'(v, 0) are finite, and 

/•OO 

/ \F(v,x)\dx < oo, (8.3) 

Jo 
/•oo 

/ li 71 ^^, x)|cte < OO 

Jo 

d.) F(u,x) is analytic in z/, for veS(T ), and fixed x > 0. 
Proof: 

In the Appendix, we prove that F(u, x) can be written as, 

I roo 



F(u, x) = - r D(v, a)e~ ax da, (8.4) 

71 Jl 



where 



1 r V —V , 



= 7r(^-l/4Ma)at[af+qT 
v y MH(i/,ia) v 7 

This result is obtained by deforming the contour in Eq. (8.2) and using the original 
representation (3.10) for M^(u,k). We note that S(u,k) is analytic for Imk > 0, 
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except on the cut k = ir; 1 < r < oo. Next we note that 

\M('\v-ia)\ > r ( Vl 2 t 2 ai)) | " \9\\Mt\ia)+R[-\v,ia)\. (8.6) 
(1 + ay 

This follows from Eqs. (6.4) and (6.13). Since (2 + a\) > l,we get 

l MH (-,^|>i + l[^±^]-|; a>l; (8.7) 

where the last term comes from lemma (6.1). Hence we have 

\M(-\v,ia)\>± (8.8) 

Finally we obtain from (8.5), with \Rev\ < 1/2, 

\D(u,a) | < 2tt(z/ 2 - ^)e-™a5, a > 1. (8.9) 

This bound guarantees the absolute and uniform convergence of the Laplace transform 
in Eq. (8.4) for all x G [0, oo), and hence all the assertions a.), b.), and c.) of our 
lemma are true. Finally, d.) is also true, given the lower bound in Eq. (8.8) and the 
uniform bound on (a u ^ 2 + oT v l 2 ) for veS{T ). 

The Marchenko Equation can now be defined as 

[•OO 

A(v;x,y) = F(v;x + y) + / duA(v; x,u)F(is;u + y) (8.10) 

with 

A(v;x,y) = 0, y < x. (8.11) 
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The integral equation (8.10) is of Fredholm type and the Hilbert-Schmidt norm 
of F is finite, 

/"CO rOO 

\F\ 2 = du dv | F{u;u + v) | 2 < oo (8.12) 

J X J X 

This follows from lemma 8.1. 

We first prove the following lemma: 
Lemma 8.2: 

For all ueS{T ) and x > 0, we have 

\F(v ] x)-F (x)\<j I ^e^\ (8.13) 
where the constant C is bounded 

C < 2 x 10 3 (8.14) 

Proof: 

From Eqs. (6.4) and (6.5) we have 

M (±) (z/, k) - M^\k) = gM[ ±] (k) + gRt\u; k), (8.15) 
where g = [y 1 — 1/4)" 1 . From the definitions of F(u,x) and F (x) we get 



1 f JU MM(u,k) M^(k), Jkx 



F(v,x) - F (x) = — / dk[ — t , ° T ^]e lkx . (8.16) 
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This gives 



F-F C = (AF) 1 + (AF) 2 , (8.17) 



with 



and 



{AF)l = ^J L e mFmh • (8 - 18) 

(AF) 2 . — / ^g^fj^l. (8.19) 

Again here L is the line k — (A + |), and — oo < A < +oo. 

In order to get bounds on (AF)i )2 , we need to separate out the terms which are 
only conditionally convergent, i.e. 0(1/ k) as k — > oo, from those which are absolutely 
convergent and hence easier to handle. 

From Eq. (5.36) and Eq. (6.5), we write 



M[ ± \k) = ^^ ) +M[ ± \k)- (8.20) 



and 



R^\u,k) = [°° da ® 1 , ( 7 N (o^ 2 + a-"' 2 ) + R?\v, k); (8.21) 
Ji [Qt±ik) 



(a ± ik) 



where 



r(± ).. A 6 



M ''W-g(T±W ; (8 - 22) 
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1=2 



Both Mi and Ri are 0(l/k 2 ) as |fc| — > oo. 
For (AF)i, we can write 



(8.23) 



(AF) 1 = (AF) 11 + (AF) 12 , 



(8.24) 



with 



+ (i-^) u m(-)m(-) 



and 



(AF)ia = — / ^A;e lfcz 



Mt ] \k)M[ +) \k) - M^\k)M\->{k) 



M (_) (A;)M(-)(z/, fc) 



(8.25) 



.26) 



The integral in (8.25) is conditionally convergent, \Mf{k)\ — > 1 as — > oo, and 
|M^(i/, fc)| is bounded from below for all k with 7mA; > 1/4. Also |Af( - )(i/, fc)| -> 1 
as \k\ — > oo. 

To obtain a bound on (AF) n we first note the following: 



dke 



ikx 



■ Mt\k) 1 
(1 -ik) ' Mt\k)M(~)(u, k) 



0. 



(8.27) 



This follows from Jordan's lemma. The integrand in (8.27) is analytic for link > 0, 
and the bracketed term is 0(1/ k) for large \k\. 
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Adding twice the i.h.s of Eq. (8.27) to Eq. (8.25), one obtains 



(AF)n ! 4~ I dh 

271 JL 



ikx 



2Mt\k) (M^(fc)-M^(fc)) 



(+)< 



rH, 



(1 + k 2 ) 



(1-ik) 



Mt\k)M(-)(u, k) 

(8.28) 

To obtain an upper bound on |(AF) n |, we first need lower bounds on Mff\k) 
for k G L. 

From Eq. (5.13) we have, 



\Mt\k)\ > 1 













a 2 




11 


— ik 




1 


- ik 


2 



(8.29) 



with |ai| = 0.654, and | a-2 1 = 0.346, and k = A + i/4, — oo < A < +oo, we finally have 

|Mo _) (A;)| > 0.255 > 1/4 ; k G L. (8.30) 

Lemma 6.1 will now give us a lower bound for M^~\v, k),v G S and k = A + if A. 
Since |AfH - M (_) | < C 2 /T 2 we get 



|M (_) (i/,A;)| > 1/4, 

for k = A + V 4 - 

Finally, we need to bound Mq +) - Mq~\ where from Eq. (5.13) 

—1ika\ 4ika,2 



(8.31) 



(1 + P) (1 + P) 



2 ' 



(8.32) 



53 



But for k = A + i/4:, \k/(l + ik)\ < 1, we obtain 



\Mt\k)-Mt\k)\< 



I -iky 
with 

64 

5 = 2 aJ + — a 2 = 2.78. 
15 

The above bounds lead us immediately to 

2MNe^ r+oo f 2 45 



|AF|ll -^^loo "Mll + (A + V4) 2 | + |l-^ 



< ° u c~ x/4 , 
\Imv\ 2 



where 

Cll = iM (2+4i) .r^< 290 . 

The bound for (AF)i 2 is easier to calculate, from (8.26) 

I „ I _ — rr/4 <.+oo if A/T^^ 

|(AF) 12 |<^- / rfA 4|M 1 (+) (A;)|+4|M 1 ( - ) (A;)|.|^' 



2vr 7-oo I 1 v yi 



For A; = A + i/A, a simple calculation gives 

M + \k) 

Mt\k) 

Using Eq. (8.22) we get 

l(AF) 12 | < 
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< 5/3. 



\Imu\ 2 



with 



The values of 6^ are given in Eq. (5.40). 

To estimate (AF) 2 we again split R^\v, k) into two terms as in Eq. (8.21). We 
write 

(AF) 2 = (AF) 2 i + (AF) 22 . (8.42) 



where 



(AF) = — [ dkc ik* [ M t\k)ri +) -Mf\k)rj- ) ] m 



and 



M^>(k)M(-)(u, k) 



(AF) 22 = — / ^[ifi^ffl; ( 8 .44) 
1 ^ 27t7l 1 Mt\k)M(-)(is,k) 



with i?[ ±) given in Eq. (8.23) and 

(a ± ik) 



r^\u,k) = dap-^-[a»/ 2 + a""/ 2 ]. (8.45) 



Equation (8.43) can be rewritten as 

(AF) 21 --J L dke I m _ ; }. (8.46) 



Here we have again used Jordan's lemma, which implies 

ikx\ IV1 r l I _ 



dke lkx [ " 1 1=0. (8.47) 
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From Eqs. (8.45) and (8.46) we get 

(8.48) 



where k = A + i/4. Using the bound (8.33) we obtain 



with 



l(AF) 21 | < (8.49) 



C2i = ^- /°°rfa(2a 1/4 )|$i 2) (a)| (8.50) 

+^ 8a 165 

-oo \a 2 - J) + A 2 + |a + 1/4 - iA||f - iA| J 



This leads to 

C2i < /° rfa« 1/4 |$f ) (a)|[^ + ^5],< (45). /°° tW /4 |$? } (a)|. (8.51) 
Ji 15 5 Ji 

Using the definition of & f\a) in Eq. (5.26), table I, and the bounds on ^\a) given 
in Eq. (6.10), one can easily get a rough numerical bound on the above integral, 
J™ daaV 4 \& 2 \a)\ < 10, and hence 

C21 < 450. (8.52) 
(AF) 22 . From Eqs. (8.44), (8.30), (8.31), and (8.23), we get 

(8.53) 
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with k = A+i/4. Using the fact that \a±ik\ > [(a^l/4) 2 + A 2 ] 1/2 , and |M (+) /M,5 } | < 
(5/3) for k G L, we have 



with 

32 & r°° . 1 

€=2 

where 



^<^Y,l daa 1 ^?\a)\(—^y-\fo, (8.55) 



/ + 0O ff?/ 
- (TT^- (8 - 56) 



(1+W 2 )^/ 2 ' 

For ct > 1, |(a/(a - 1/4) | < 4/3 and hence 



oo 5 

c 2 2<^J2 A ^Um £ (8-57) 



37F £=2 



with 



A e = J™ daa 5/4 - e \$¥\a)\ (8.58) 
A simple numerical estimate will give 

c 22 < 10 3 . (8.59) 

This completes the proof of our Lemma with C = c n + c 12 + c 2 i + c 22 < 2 x 10 3 . 

Lemma (8.2) guarantees that for v G S, \F\ = \F Q \ + O(-^j). Indeed we have for 
the Hilbert-Schmidt norms 

\F\ < \F \ + \F- F \ (8.60) 
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However, from the lemma 

poo poo 

iv\F(u,u + v) -F (u + v)\ 2 (8.61) 



poo 


du I 


J X 


j . 


c 2 


roo 


|z/| 4 


Jx 


c 2 


roo 


M 4 


Jo 



X 



< T-TT I du [ dve 2 



< i—rr I we 2 dw. 
Thus 

2C 2 

|^-^o|<7^<^- (8-62) 
We have the exact expression for F (x) given in Eq. (7.5) and we can calculate 
|F | exactly. The Hilbert-Schmidt norm for |F | is 

roo roo 

\F \ 2 X = du dv[F (u + v)} 2 . (8.63) 

Jx Jx 

Note that here x appears as a parameter (see Eq. (8.10)). As x increases the norm 
of Fq tends to zero. From Eq. (7.5) we get 

POO 

\F \ 2 = / dw[\ 2 w + 2 lambda Q \ lW 2 + Ai^le -2 "', (8.64) 
Jo 

with A and Ai given by Eqs. (7.6) and (7.7). After some algebra we obtain, 

\Fo\l =0 = l + ^ , (8.65) 

with A x = —0.590. Thus the Hilbert-Schmidt norm for x = is slightly bigger than 
one, however, it is easy to show that for x > x Q = 0.1, 

\F \ < 1. (8.66) 
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Hence the iterative series for both F and F D converges for a; > 0.1. But to prove 
the existence of A(u; x, y) for all x > 0, we will use Eq. (8.61) and proceed in another 
way. 

In section VII we gave an explicit solution of the integral equation, 

A = F + A F . (8.67) 



We can write 



I = (l-F )'\ (8.68) 



and 



A — (1 — F Q ) F Q — I F a . (8.69) 



This leads to 



A + 1 = I . (8.70) 

The kernel A a (x, y) is given explicitly in Eq. (7.9) and (7.18-19), and thus I (x,y) is 
known. 

The full Marchenko equation (8.10) can now be written as 

A = F + AF. (8.71) 

We define / as 

/ = (l-F)-\ (8.72) 
59 



and 

A=(1-F)- 1 F. (8.73) 
We prove that both I and A exist and in (8.70) have a finite norm. We have 

A+l = I = (1 - F)-\ (8.74) 

To show that I (and A) exist, we note first, 

/ = (1-F)- 1 = (1-F - A)" 1 , (8.75) 

where 

A = F — F Q . (8.76) 
Using the fact that (1 — F )I = 1 we get 

/ = [(1 — F )I {(1 — F ) — A}] -1 , (8.77) 
= (l-IcA)- 1 ^. 

Using Eq. (8.74) we get 

A = (l-I A)- l I -l (8.78) 
= (1-I A)-\A + I A). 



where A Q + 1 — I a . 



60 



Next we define the operator K as 

K = I Q A = (A + 1)A, (8.79) 

and obtain 

A = A + (1-K)- 1 .(1 + A )K. (8.80) 

The Hilbert-Schmidt norm of A is small. Indeed using Eq. (8.61) we have 

|A| 2 = |F-F | 2 <^< (8.81) 

However, K = A a A + A, and we get 

\K\ < \A \.\A\ + |A|. (8.82) 

A Q is known and \A Q \ 2 < 5, thus 

M < « 1, (8.83) 

for all z/ G S(T ). Thus the inverse(l — if)" 1 is given by an absolutely convergent 
series, 

oo 

(1- K)- 1 = Y,K n , (8.84) 

n=0 

and has a bounded norm, |(1 — K) _1 \ < 2. 
The fina result for A is, 

A = A + H + A H, (8.85) 
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where 

oo 

H=Y,K n ; (8.86) 

n=l 

and \H\ < 2\K\ « 1. Indeed we have \A - A \ < C = j\v\ 2 . 
The kernels (z/; x + y) can be written as 

[•OO 

K(u; x, y) = F(u; x, y) - F (x + y)+ duA (x, u) [F(u; u + y)- F (u + y)} (8.87) 

Jx 

The properties of the kernel K{y \ x, y) are similar to those of F(v; x, y). We have 
the following: 
Lemma 8.3: 

K(u;x,y) is for y > x > 0, a.) analytic for v G S(T Q ); b.) different iable in 
both x and y; c.) analytic for Rex > and Rey > when v e S{T ); and d.) 
\K(v;x,y)\<^e-(-^. 

Proof: 

These results follow from Eq. (8.88), lemma (8.1), and the exact result (7.11) 
for A (x,u). We note that the denominators appearing in the expressions (7.18) and 
(7.19) for B(x) and C(x) do not vanish for Rex > 0. Both B and C are thus analytic 
in the half plane Rex > 0. 

The full expression for A{y\ x, y) is 

[■CO 

A(y; x, y) = A a (x, y) + H(y; x, y) + / A (x,u)H(is;u,y). (8.88) 

Jx 
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where 

oo 

H(v;x,y) = Y,KW(v,x,y), (8.89) 

n=l 

and 

K {n \v]x,y) = f dui... I du n - 1 K(u;x,u 1 )K(u;u 1 ,U2)....K(u;u n - 1 ,y) (8.90) 

J X J X 

The series in (8.89) is absolutely and uniformly convergent for y > x > 0, and all 

v e S(T ). 

The properties of A{y\ x, y) can be summarized in the following lemma: 
Lemma 8.4: 

a.) For y > x > 0, A{y\x, y) is analytic in v for v G S(T ); b.) A{y\x, y) 
is differentiable in both x and y, y > x. Also A(z/, 0,0), and [-^(A(is; x, x)] x=0 are 
finite; c.) For fixed z/, v G <S(T ), A(z/;x, y) is analytic in a; and y for Rex > 0, 
-Rey > Rex > 0; and d.) For all v G S(T Q ), we have the bound 

\A(u; x, y) - A Q (x, y) \ < ^.e~^\ (8.91) 

\ v \ 

Proof: 

These results follow immediately from Eq. (8.89) and lemmas (8.1) and (8.3). 
The bound (8.91) follows from the bound (8.13) of lemma (8.2). The constant C is 
certainly such that, C < 10 4 , which is sufficient for our purposes at this stage, but 
can be improved with more careful estimates. 
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The next step is to define two functions U(u;x) and /^(v; k,x) as follows: 

U(v;x) = -2—A(i>:x,x),x>0; (8.92) 
ax 

and 

f ± (is;k,x) = e* ikx + dyA{v-x,y)e^ ky . (8.93) 

Jx 

Without recourse to the standard methods of inverse scattering, one can prove 
directly the next lemma. 
Lemma 8.5 : 

For any v G S(T ), f ± satisfy a Schrodinger equation with U (u, x) as the potential 

-^ + U(v;x)f ± = k 2 f ± . (8.94) 

Proof: 

From a.) and b.) in Lemma (8.4) it follows that U(u;x) is analytic in v for 
v E S(T ) and x > 0. Similarly from Eq. (8.94) it follows that f^{v;k,x) with 
x > 0) and Imk > is also analytic in v in the truncated strip. Similarly with 
Imk < is analytic. The same is true for -^(f^), and -^(f^) since absolute and 
uniform conveergence allows us to differentiate under the integral sign in (8.93). 

In the next section we will prove the validity of Eq. (8.94) on the line v = it, 
t > T Q . Hence by analytic continuation, the Schrodinger equation (8.94) holds for all 
y e S(T ). 
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In the Appendix will give a more direct proof of Eq. (8.94) and also show that 
one does indeed recover the original Jost function from the potential U(u; x). 
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IX. The Case, v = it. 

For purely imaginary u, our S-matrix, S(v,k), is unitary and satisfies all the 
properties needed for the standard inverse scattering methods of Gelfand, Levitan 
and Marchenko to be applicable. We sketch some relevant results in this section. 

First, we define S as: 

«>•*• < 9l > 

For real k, it follows from Eq. (3.10) that 

[M i+ \it, k)}* = M (+) {it, -k) = M ( ~\it, k); (9.2) 

S(it, k) satisfies all the conditions given in Faddeev's 14 review paper, which are 
sufficient to guarantee that the Marchenko equation will lead to a unique real potential 
U(it,x). We can easily check that for real k, 

\S(it, k)\ = S(it, 0) = S(it, oo) = 1; (9.3) 

and 

[S{it,k)]* = S{it,-k). (9.4) 

The number of discrete eigenvalues of S for fixed t > n 2 is at most one (see Section 
VI). In the physicist's language, we have either one bound state or one antibound 
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state. This is evident from Eq. (6.28) and the fact that £(it + 1/2) is real. We will 
discuss this point in more detail at the end of this section. 
The Marchenko kernel is now given by 

I roo 

F(it, x) = — dk[S(it, k) - l]e tkx ; x > 0; (9.5) 

for the case, 

t(it+l)>0. (9.6) 

and 

F(it, x) = - dk[S(it, k) - l]e ikx + c e~ T ° x . (9.7) 

7T J-oc 

for the case where £(it + 1/2) < 0, i.e. with a bound state at E — —r 2 . Both Fourier 
transforms in Eqs. (9.5) and (9.6) are convergent in the mean since [S — 1] = 0(1 /k) 
for large k. Also it is clear that F(it, x) is real. 

As noted previously £(it + 1/2) = 0(t p e~^) and hence small for, t > w 2 , this 
makes r G << 1/4 and the bound state is very shallow for t > n 2 . One can now move 
the contour of integration up in both Eqs. (9.5) and (9.6) to obtain 

F(it, x) = / dk[S(it, k) - l]e ikx ; x > 0; (9.8) 

for both cases. Here L is the line Irak = 1/4. In the case of Eq. (9.6), the contribution 
from the pole at k — ir Q exactly cancels the second term on the r.h.s., see ref. 22. 
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The solution of the Marchenko equation, A(it;x, y) exists, is real, and differen- 
tiable for y > x > 0. The resulting potential, U(it,x), is real, continuous for all 
x > 0, and 0(e~ 2x ) for large x. 

We close this section by calculating the position of the bound state or antibound 
state for fixed v — it, t > T Q . 

Rewriting Eq. (6.26) for v = it, we get 

M { ~\it,k) = 2i(it+-) + M^\k)+gM[~\k) (9.9) 
+ ^fda (^.{j^y - ^} *f»(a)[co S (i loga)]) , 

where 

3 = WTTJi (9 - 10) 

Now, £(it + 1/2) is real and exponentially small for large t. The three other terms on 
the r.h.s. of Eq. (9.8) are all 0(k) for small k. 
From Eq. (6.13) we have, 

M^\k) = -tk(2 + ai ) + 0(k 2 ), (9.11) 

with 

Gl = -1 + 8^(1) = -0.6543. (9.12) 
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Thus for t > T , the one zero of \it, k) will occur only when k = ir and 

where the 0(l/t 2 ) term is real, and the same fore the 0(r 2 ) term. Note that M^~\it) 
and m[ \ir) are both real as is the integral in (9.8) for k = ir. 
The resulting potential, or one parameter family of potentials, 

U(it,x) = V(g,x), (9.14) 

has a remarkable property as t increases, t > T D . It will have exactly one bound state 
when £(it + 1/2) < 0, with energy E Q = — r 2 , 

^ = l(2?")Vo(f)P +0 «« ft+1/2 )'' ) - <9 ' 15) 
Then as we pass a Riemann zero and £(it + 1/2) > 0, there will be no bound state 
until t reaches the next Riemann zero. 

As t — > +oo, the potential, U(it,x), presents us with a seemingly puzzling situa- 
tion. The bound state, i.e. a point spectrum of one, appears and then as t increases 
disappears, with this process repeating as t increases until as t — > oo we reach V (x) 
which has no point spectrum. 

Schwinger's theorem relating the number of bound states to the number of nodes 
of the zero energy regular solution, <f>(it; 0; x) , defined in Eq. (2.8), is instructive for 
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the present case. 

From ref. 13, we have an integral equation for 0, 

sinkx f x sinkix — x') T r , .. , . , , . 

<p(it; k; x) = — — + / K - '-VMWit; k; x')dx', (9.16 

k Jo k 

where clearly (j>(it; k; 0) = 0. The zero energy <p is given by 

(f)(it; 0; x) = [1 + / V((x')(f)(it; 0; x')dx']x - / V{x')<f){it;0,x')dx' (9.17) 

Jo Jo 

For x not large, x < T a , we can approximate by the t — > oo solution 

0(^;M) = O (M) + O(^), (9.18) 

where O is defined as in (2.8) but with J 1 * 1 replaced by and M ± replaced by Mf. 
It is easy to check that O (O, x) is positive for x not large, and O (O; 0) = 0. But from 
Eq. (9.17) we see that the large x behavior is 

<p(it;0;x) -> [C(t) + o(l)]+ x[l + / V(ar')0(i*; 0; x')^ + o(l)] (9.19) 

Jo 

where 

C(t) = - x'V(x')(t)(it;0;x')dx'. (9.20) 
But under the integral sign we can replace V by V Q and (f> by O , and obtain 

/oo 1 

xV o (x')0 o ^; 0, x')dx' + O(-). (9.21) 
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where V and O are known exactly from section VII. It is a simple matter to check 
that 

C(t) > 0, t> T Q . (9.22) 
Next in ref. 13, we have the result relating the Jost function to 0, and it gives 

M { -\it; k) = 1 + j°° dx'Viit; x')(f)(it; k; x'). (9.23) 

Jo 

Taking the k — > limit we have from Eq. (9.19) 

<j>(it;0;x)^C(t) + 2£(it + -)x, x -> oo. (9.24) 

We can only have a node in if + |) < 0, otherwise there is no node and no 
bound state. For t > T Q and large, the node occurs at large values of x, 

This discussion shows that while our asymptotic estimates for /, 0, and V, are 
good for low values of x, x < T a , one cannot use them for large values of x except in 
estimating integrals as in Eq. (9.21). 

The well established results for a unitary S and real potentials now guarantee that 
Schrodinger's equation holds for f ± it; k;x) and U(it;x), i.e. 

-^f ± (tt;k;x) + U(tt;x)f ± = k 2 f ± , t > T . (9.26) 
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We are also guaranteed that 

f ± (it; k; 0) = M^\it\ k), (9.27) 
where is the original Jost-function we started with 

MW(it; k) = 1 - (t 2 + 1/4) f ^) a^a^ + ^] ^ 

Thus by analytic continuation the Jost solutions f ± (u;k;x) given in the previous 
sections will also give the original Jost functiona, i.e. (9.28) with it replaced by u, 
and v G S(T ). 
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X. Asymptotic Expansion in Powers of g 

In this section we carry out the asymptotic expansion of the kernels F and A and 
the potential V in powers of g, where g = (v 2 — 1/4) _1 , and ueS(T ) with T Q > 10 3 . 
We start with the definition of the Marchenko kernel F(u; x) given in Eq. (8.2), 

F{y;x) = ^-l dk[S(u,k)-l}e ikx ,x>0, (10.1) 

where L is the line Irak = 1/4. 

Using the asymptotic expansion of M^(u, k) given in Eq. (5.47), we have 

N 

J2g n Mi + \k)+g N R^(u,k) 
-1 = 2^ 1, (10.2) 

n=0 

where i?^(|nw, fc) are both 0(g) for small g. 

Next we recall the lower bound obtained in section VIII for M^~\k) when k = 
A + 2(1/4), and —00 < A < +00. This is given in Eq. (8.30) 

I M(~\k) |> 0.255, *eL. (10.3) 

From lemma 6.1, we get also a lower bound on \M^~\v, k) \ for ueS(T ), and /ceL given 
inEq. (8.31) 

I M { - ] (u,k) |> 1/4 , T > 10 3 . (10.4) 
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The last two bounds guarantee that the denominator in Eq. (10.2) does not vanish 
for any keL and ueS(T ). We can then proceed to expand [S(u, k) — 1] in powers of 
g for any keL and get 

N 

[5(1/, k) - 1] = £ g n H n {k) + 0*4%, *) (10-5) 

n=0 

where from Eq. (10.2) we get 



tf o(W )-l = «-l ; 



;io.6) 



M^\k)M[-\k) 
Mo~~\k) V±V±l VV Mo~\k) 



H i(k) = TttTzi^W " ~ ~ >, (10-7) 



and 



S * {M<+>(fc) - ^ - } W (10.8) 

Mo (k) Mt\k) Mt\k) 

MM(k)[Mt\k)}\ 

[Mo~\k)\ 2 



with similar expressions for H n (k), n > 2, which we will not need in this paper. The 
remainder term is O(g) as g — > 0. All H n (k) are rational functions of fc. 

Eq. (10.5) immediately gives us the asymptotic expansion for the kernel F(z/; x), 

F(i/; x) = F (x) + gF 1 (x) + g 2 F 2 (x) + ... + g N F { R N) (g; x) (10.9) 
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where 

F n ( x ) = f dkH n {k)e ikx ,x> 0. (10.10) 

This last Fourier transform is conditionally convergent since \H n (k)\ = 0(1 /k) for 
large k (note that L is the line Irak = 1/4). 

Next we stress that all the H n (k) are analytic for Irak > 0. The denominators in 
Eqs. (10.6)-(10.8), do not vanish in Irak > 0, except at k — 0. 

This also holds for H^\g, k). Thus we can shift back the contour in Eq. (10.10) 
to the real /c-axis and obtain 

Fn(x) = — / dkH n (k)e lkx ,x> 0. (10.11) 

for real k, it follows from Eq. (5.49) that [M^(k)\* = M^(k), and that M^(-k) = 
M^-\k). This leads us to, 

H*(k) = H n (-k), fork real. (10.12) 

Thus it immediately follows from Eq. (10.11) that all F n (x), n — 0, 1, 2, N, are 
real functions. However, F^\g,x), is certainly not real for veS(T ). 
We have explicitly calculated F (x) in section VII, and obtained 

F (x) = \ e~ x + \ x xe~ x , (10.13) 

with A G and Ai real and given in Eqs. (7.6) and (7.7). 
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One can also easily calculate explicitly Fi(x) by contour integration. 

Fl ( x) = -1 rdk\*^W - M M k ) M t\k) ** (10 14) 

2nJ-oo [ M^'\k) [Mt\k)] 2 

The result is 

F 1 {x) = {Y,<r n x n )e- x , (10.15) 

n=0 

where the constants a n are explicitly given as functions of ai, and 6j, j = 1, 4. Here 
we will only give the numerical value of the a' n s 

a = 26.5228, (10.16) 

cti = 1.7901, 

a 2 = -9.3291, 

(7 3 = -2.3582. 

The result for F 2 (x) will be similar, 
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F 2 (x) = e- x Q2P n x n ). (10.17) 



n=0 



We will not give its explicit value as it is not needed. 

The Marchenko Equation (8.10), for ueS(T ) can be written in operator form 

A = F + AF (10.18) 



76 



Writing 

A = A + gA 1 + g 2 A 2 + g 2 A%\ (10.19) 
and using the expansion for F given in Eq. (10.9) we get by comparing terms, 

A = F + A F , (10.20) 

an equation we solved explicitly in section VII. In addition we have 

A 1 = (Fx + A F 1 ) + A 1 F , (10.21) 

as the integral equation for A ± , and 

A 2 = (F 2 + AiFi + A F 2 ) + A 2 F , (10.22) 

for A 2 . It is obvious that the integral equations for A n ,n = 0,1,2,... all have the 
same kernel F . Thus they are all explicitly solvable. Given our solution, A Q , for Eq. 
(10.20) we get, 

A Q = (1 - F )- l F 01 (10.23) 

and hence 

(1-Fo)- 1 = A + 1. (10.24) 
This leads to solutions for A 1: A 2 , etc. with 

A 1 = F 1 + 2A F 1 + A (A F 1 ), (10.25) 
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and, given Ai, we can now get A 2 explicitly as: 

A 2 = F 2 + A 1 F 1 + 2A F 2 + A (A 1 F 1 ) + A {A F 2 ). (10.26) 

It is now obvious that all the A n 's are real and continuously differentiable, for 
y > x > 0, since from Eq. (7.11), 

A (x,y) = [B(x) + (y- x)C(x)]e^*\ (10.27) 

with B(x) and C(x) given by Eqs. (7.18) and (7.19) and B, C are 0(e~ 2x ) as x — > oo. 
The kernels F n (x) are of the form 

2n+l 

^n(x) = e^[5: ofV]. (10.28) 

The potential C/ {y \ x) is given by 

U(u; x) = -2^-A(u; x,x),x>0. (10.29) 
ax 

Using the variable g = (y 2 — 1/4) -1 , we write 

C/"(i/; x) = V(g;x). (10.30) 
The expansion of A in powers of g, given in Eq. (10.19), now gives us, 

V(g; x) = V (x) + gV^x) + g 2 V 2 (x) + ... + /V^, x), (10.31) 
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where 



V n (x) = -2— A n (v;x,x), 



(10.32) 



V™tg,x) = -2—A<P(v;x,x). (10.33) 

and V^ N) is 0(g). 

All the V^j's are real, continuous for x G [0, oo), and 0(er 2x ) for large x. V^ N \g, x) 
is complex for veS(T ) but also continuous and 0(e~ 2x ) for large x. 

It is now clear why one can refer to " g" as a coupling constant specially for large 
values of Imu, i.e., small g, g — 0(\lmv\~ 2 ). 

We have already calculated the first term in the expansion, V (x), and it is given 
explicitly in Eq. (7.21). Later we will need to have V\(x) and we proceed to calculate 
it here. 

From Eq. (10.25) we get, 

POO 

Ai(x,x) = F 1 (2x) + 2 dzA (x,z)F 1 ( y z + x) (10.34) 

JX 

POO POO 

+ / dzi dz 2 A (x;z 1 )A (x;z 2 )F 1 (z 1 + z 2 ), 

Jx Jx 

where A Q and Fi are given in Eqs. (10.27) and (10.25). It is now evident that Ai(x, x) 
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is continuously differentiable and 



2^1 



dx 



(x,x), 



(10.35) 



can be explicitly calculated in terms of B(x), C(x), and the real constants aj. The 
resulting V\{x) is continuous, finite at x — 0, and V\ = 0(x 3 e~ 2x ) for large x. We do 
not write down the full result, but exhibit a graph of V\{x) in figure 1. 

In closing we comment on the asymptotic expansion (10.31) for V(g;x). It of 
course is not convergent no matter how small \g\ is. This follows from the fact that 
the expansions for M ± (u; k) are also divergent. There is an essential singularity at 
g — 0. The constants xT\ ^ = 1? ■■■■> 2n + 2, given in eqs. (5.32) and (5.33) grow fast. 
However, Eq. (10.31) can still give an extremely good estimate for V(g; x) as long as 
N is 0(1). Indeed it is possible to get a uniform bound on V^ N \g,x) which is 



where C(N) grows fast with N. For the purposes of this paper we need at most N = 2 
or 3. For T = 10 4 , (C(2)/T D 2 ) « 1. This makes {g 2 V^) smaller than (rA)e- x , 
and thus V + gV\ + g 2 V 2 , give an excellent estimate of V(g, x) for all veS(T ), and x 
not large. However, this estimate is not good for large x where both V and the error 




T I9 C ' 



Imv\ 2 



x > 0, 



(10.36) 



are small. 
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Finally, we will need an important result on the phase of V R near the critical 
line v = it. 

As we have shown in section IX, F(it;x), A(it; x; y) , and U(it;x) are all real, for 
xe[0,oo) and y > x. In addition, in the asymptotic expansion given in Eq. (10.31) 
all the coefficients V n (x) are real. But the remainder term, V R f \g,x), is in general 
complex for v = uj + it, and w^O. However, for uj = 0, we again have reality, 



In Section VIII we proved that both A{y\x, y) and U(v;x) are analytic in u for 
v G S{T ). Hence so is U^\v\x). This leads us to the following lemma. 
Lemma 10.1 : 

For ueS(T ), and v — 8 + it, we have 



u^( i t-,x) = vl»\-(t> + 1 -r-,x) = (v^y. 



(10.37) 



\Im(Vjf\g;x))\ = \ImU { R N \5 + i t-,x)\<C(x).\5\+0(S 2 ), 



(10.38) 



where 



C(x) < -ie- 



X' 



(10.39) 



Proof: 



[ImU R (u;x)] = Im 



dU R "\v;x) 



(10.40) 



The derivative is finite, and setting v' — 5 + it, 5 « \, we get equation (2.38). 
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The bound on C(x) follows from our previous estimates. We will not give the 
proof here. 
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XI. The Zeros of M' \v\ k) for Fixed k 

In this section we shall study the properties of the infinite set of zeros, v n (k), of 
M(~\v, k) for fixed k. We prove three lemmas for {u n (k)}. 

For convenience and without loss of generality, we set k = ir, and r > 0. We 
write 

M(v-T) = M { ~\v,ir), (11.1) 

and 

M n {r) = M { -\ir). (11.2) 

It is clear from the equation defining M(u, r), i.e. Eq. (3.10) with fc = ir, that 
M(v,t), fixed r > 0, is an entire function of v with order the same as £(f + 1/2), 
i.e. order 1. Hence M(u, r) has an infinite set of zeros, ^(t), with (^(r)! — > oo as 
n — > oo. 

Next it follows from section VI that, ^„(r), will all be outside the truncated critical 
strip, S(T a ), for r > T^e" 2 ^. As we decrease r, i^ n (r), will start appearing in S(T ) 
for r < 0(T p e-?f). 

Our first lemma is the following: 

Lemma 11.1: 
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As r -> 0, 

lim T ^ u n (r) = v n , (11.3) 
where (z/ n + |) is a zero of the zeta function i.e. 

£K + 1/2) = 0. (11.4) 

Proof: 

From the asymptotic expansion of M{y\r) given in Eq. (5.47), we have 

M{v- r) = 2£(v + 1) + M (r) + gM^r) + £ 2 M 2 (r) + g 2 R 2 (u, r), (11.5) 

where 

/•OO ? f 1 1 ^ 

R 2 (v,t) = 1 da(^?\a) [j^^ - -^)[a% + a^}. (11.6) 

In getting Eqs. (11.4) and (11.5) we have used Eq. (5.46) for the ^-function. 
By definition we have, 

M(u n (r),r) = 0. (11.7) 

Hence we get 

-2£K(r) + i) = M^+^Mxto + i^MaO-) (11.8) 
+ ^(t)£ 2 K(t);t), 

with 

9 - (T) = i*hr (1L9) 
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Now all the terms on the r.h.s. of Eq. (11.8) are 0{t) as r — > 0. Hence we get 

lim T ^(u n (r) + 1) = £K(0) + ^) = 0, (11.10) 

and therefore 

^n(O) = !/„. (11.11) 

The next lemma gives us as estimate of {vu{t) — ^n(O)) as r — > 
Lemma 11.2 : 

If v n is a first order zero of £(f + 1/2), then as r — > 

K(r)-i/ n ) =0(r), (11.11) 

and if v n is of order p then 

K(r)-z/ n )=0(r^). (11.12) 

Proof: 

Since £(f + |) is entire we can write for any z/ n 

£K(r) + \) = £K + \) + ^l = „ n (v n (r) - v n ) + 0\{v n {r) - v n )\ (11.13) 

But the first term on the right is zero, and (£') v=Vn ^ for a first order zero, and we 
get from Eq. (11.7) 

- 2(0^>n(r) - i/„) = r[(2 + Ql ) + 0(<7„(r))] + 0(r 2 ). (11.14) 
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This gives our proof for a first order zero. 

For a zero of multiplicity p, we have by definition (^) u=Un = 0, for j = 1, 2, ...,p— 
1, and (^) u=Un 7^ 0. Hence we get, 

- 2[^X =Vn {v n {T) - vf = r[(2 + fll ) + 0(</ B (r))] + 0(t 2 ). (11.15) 
This leads to 

(u n (r)-u n )=O(T 1 / p ),r^0. (11.16) 

So far we have shown that every z/„(r) approaches a Riemann zero as r — > 0, but 
have not established the converse, i.e. that any u n is the limit of a u n {r) as r — > 0. 
To do this we first define a rectangular region i?(T , T) as follows: 

i?(T , T) = {v\ - 3/2 < Rev < 3/2; T < Imv < T}; (H-17) 

with T»T Q . 

We now prove our third lemma. 
Lemma 11.3 : 

Let A^(T G , T) be the number of zeros of \) for veR(T , T), and N M (T , T; r) 
be the number of zeros, ^ n (r), of M(z/, r), with u n (T)eR(T ,T), then for sufficiently 
small r 

|iV M (T , T; r) - iV 5 (T , T)| < 1. (11.18) 
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There exist a small interval in r, < r < r (T), such that 

N m (T ,T;t)=N^T ,T). (11.19) 

Proof : 

We start with the standard expression: 

5 2nifr R X M(z/;t) £(z/ + 1/2) J ' V ; 

where T R is the boundary of the rectangle R. We also choose T, and T G such that 
they both lie between the abscissa of successive zeros v n \ i.e. Imv ni < T < Imv Nl+ i, 
and Imu rio < T Q < Imv rio+ i. Thus Yb never has a zero of £ on it. The prime in 
(11.20) denotes {d/dv). We follow the method used to prove theorem 9.3 in ref. 18. 
Using the symmetry of £ in z/ we have 

iV ? = - / 2 / ci4+- / di/l = -{Aargt (u+12 }; 11.21 

where A denotes the variation from iT Q to | + iT Q then from (| + iT ) to | + iT, and 
thence to %T . 

But from Eq. (11.5) we get 

{ M(„; T )/{(„ + !)} = !+ (2 ;; ( f + + ) CgT + OM; (11.22) 

On the horizontal parts of T^, we can use theorem 9.7 of ref. 18 to obtain a lower 
bound on |)|. Indeed, there is a constant, A, such that each interval (T, T+ 1) 
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contains a value of t for which, 

\C^+^)\>t- A ,- 3 -<u J < 3 -, (11.23) 

where v — uj + it. On the vertical parts of r B ,\((i / + |)| is obviously bounded from 
below. Using the standard asymptotic expression for r(| + |) as t — > oo, we finally 
get 

|[ar^M(i/;r) - ar^(i/ + i)]| < ^^|, (H-24) 
We can choose r = T~ A ~ N (e~¥~) with iV > 3, and obtain 

A(argM(u, r) - antf(i/ + \)) < (11.25) 

Hence as T — > oo, we get 

jV M -jV e = (11.26) 

This completes our proof. 

We stress one important fact that is a consequence of the results of this section. 
Namely, we are now not limited to the study of the zero energy zeros, z/„(0). One can 
consider the case for small enough r, but r > 0, and obviously if there is an interval 
< r < r n for which [Reu n (r)] = 0, this will be sufficient for the validity of the 
Riemann hypothesis. In the next section we will see the importance of this remark. 
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XII. The Potential V(g, x) and the Riemann Hypothesis 

Following the notation of the previous section we define the Jost solution f(g;r, x) 

as 

f(g-,T;x)=f(-\g;iT,x), (12.1) 

with k = it, and ueS(T ). We also recall the result given in Faddeev's review, ref. 
14, 

roc 

\f _ e — | < Ke -™ u \V(g,u)\du, (12.2) 

J x 

hence / = 0(e~ TX ) as x — > oo, r > 0. 

We write the Schrodinger equation for / and /*, 

-d 2 f 

dx 2 

and 



+ V(g,x)f = -r 2 f, (12.3) 



-d 2 f* 



+ V*(g,x)f* = -T 2 f*. (12.4) 



dx 2 

Multiplying the first equation by /* and the second by /, integrating from x = to 
x = oo, and subtracting the two equations we get 

roc 

2t / dx\f(g;T;x)\ 2 [ImV(g,x)] = M(iy,T)A*(u,r)-M*(u,T)A(u,T), (12.5) 

Jo 



where 



w v df(g;T;x) 

A(^,r) = [ — ] x=0 , (12.6) 
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and g = (v 2 — 1/4) 1 . The derivative exists for x — > in our present case, since V(g, 0) 
is finite. One can also check this from the expression for / in terms of A(v; x, y), 

poo 

f(g; r, x) = e~™ + / dyA(u; x, y)e~^ '. (12.7) 

J X 

From section VIII, we know that A(v;0,0) is finite, and also (dA/dx), for x > 0, 
y > x exists and is integrable. 

Next we set v = u n (r), and g = g n (r) in Eq. (12.3), and we get 

/ dx\f(g n (r);r;x)\ 2 [ImV(g n (r);x)\ =0 (12.8) 

J o 

This last integral is absolutely and uniformly convergent since V = 0(e~ 2x ) as x — > oo. 
Thus we can take the limit r — > and obtain, 

/ rfx|/(^;0;x)| 2 [/m\/(( ?n ; a ;)] = 0. (12.9) 

J o 

Now l / ((?,a:) has an asymptotic expansion in g. 

V(g; x) = V (x) + gV 1 (x) + g 2 V 2 (x) + g 2 vg\g, x). (12.10) 

Also from the expansion for A(u; x, y) we get for the Jost solution 

/(<?; r; x) = / c (r; x) + ^A(r; x) + g 2 f 2 (r; x) + (? 2 / R ((?; r; x). (12.11) 

Both Vft and are 0(g) and have bounds in x. It is sufficient, for the validity of 
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the Riemann hypothesis for a constant c > 0, c = 0(1), to exist such that, 

/ dx\f o (x)\ 2 V 1 (x) = c>0; (12.12) 

Jo 

where f Q (x) = f o (0,x). 

To prove this last statement, we first note that setting v n = u n (0) we have 

v n = LO n + it n , t n > T Q ; (12.13) 

and 

Img " = M^t+W ; w "- 1/4 ' (12 ' 14) 

The vanishing of img n with t n > T Q > 0, implies w n = and hence s n = 1/2 + it„. 
Next we write from Eq. (12.10), 

ImV{g n -x) = (Img^V^x) + (Img 2 n )V 2 (x) (12.15) 
+ (Img 2 )(ReV^(g n] x)) + {Reg 2 n ){ImVg\g n] x)). 

Also in the integrand in (12.0) we can write 

\f(g n ;0;x)\ 2 = \f(0;0;x)\ 2 + O(±). (12.16) 

In Eq. (12.15) we note that (Img n ) = 0(-&) and lmg 2 n = O($0, winLe(Reg%) = 
0(4-). On substituting Eqs. (12.15) and (12.16) in (12.9), and using the assumption 
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(12.12) we see at first that consistency requires u: n to be small, i.e. u n = O(-k-). Here 
we use the fact that \V^\ = 0(g) for small g. However, we have more information on 
V'r 1 and specifically its phase for small u n . This was given in lemma 10.1, where it 
was shown that ImV^ = 0(u) for small u> and that is real for uo = 0. Given the 
bound on ImVjp from this lemma we see that the leading contribution from (12.15) 
to Eq. (12.9) must come from [(Img n )Vi\ and cannot be cancelled by the other three 
terms. We have 

(Img n )c + (ImgDX, + (Img 2 n )X 2 + (Reg 2 n )X 3 = 0. (12.17) 

with 

|*i| < 2 / \f (x)\ 2 \V 2 (x)\dx; (12.18) 

Jo 

\X 2 \ < 2 r\f (x)\ 2 \ReV^(g n ;x)\; 

Jo 

\X 3 \ < 2 r\f (x)\ 2 \ImV^(g n ,x)\dx 

Jo 

< ^ r \fo(x)?e-*dx 

t n Jo 

Hence we have constants Bj such that \Xj\ < Bfj = 1,2, and |X 3 | < (^)B 3 . 
Note that in the last inequality we used lemma 10.1. We take T Q large enough such 
that W « c, j = 1, 2, 3; and with c = 0(1). Next we rewrite Eq. (12.17) as 

1 o 

(Img n )[c + 2(Reg n )X 1 + 2(Reg n )X 2 + (Reg 2 n ).t n X 3 ] = 0; (12.19) 
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where now \X 3 \ < 2B 3 . The term in the square bracket cannot vanish, and hence we 
obtain 

PCX) 

[Img n \c = [Img n \ / \f {x)\ 2 V 1 {x)dx = 0, (12.20) 

J o 

and thus, if c 7^ 0, then for all t n > T Q , we get 

lmg n = 0;or s n = ^ + it n . (12.21) 

We have already calculated f (x) exactly in section VII, and V±(x) in section X, 
and we can compute the integral in (12.12) directly, the result is 

/ \f o (x)\ 2 V 1 (x)=0. (12.22) 

J o 

This can be checked numerically, and indeed can be rigorously proved. So we have 
no information on (Img n ) from (12.17). However, the proof of (12.22) suggests to us 
how we can proceed further. 

To prove (12.22) we use the Schrodinger equation and the expansions (12.10) and 
(12.11). We obtain 



d 2 fo, 



dx 2 

and 



(r, x) + V (x)f (T, x) = -r 2 f (r, x), (12.23) 



^(r, x) + V (x)h(T, x) + V^Ut, x) = -r 2 /i(r, x). (12.24) 



93 



Multiplying the first equation above by fi and the second by f a , integrating from zero 
to infinity and subtracting, we have 

/ [f (r,x)] 2 V 1 (x)dx = -[/i(r,0)/^(r,0) - / (r, 0)f[(r, 0)], (12.25) 

Jo 

where the prime denotes (d/dx). 
By definition, we have 

f o (r,0) = M o (r), (12.26) 

and 

/ 1 (r,0) = M 1 (r). (12.27) 

Both M d (t) and M x (r) are O(r) for small r and vanish as r — > 0. We then get after 
taking the limit, 

PCX) 

/ |/ o (0,x)| 2 Vi(a;)tfa; = 0. (12.28) 

Jo 

The crucial factor here is the fact that M(v, 0) = 2£(z/ + 1/2), and for large \Imv\, 
veS(T ), M(z/,0) = 2^(i/+l/2) = 0(e^^). This forces all the coefficients, M n (r), 
in the asymptotic expansion of M(z/, r) in powers of g to vanish as r — > 0. The culprit 
is the factor T(| + 1/4) in the Eq. (3.4) which relates £(s) to C( s )- We will also see 
below how this fact hinders us in treating the case r^O, but r small. 

From the results of sec. XI, it is evident that it is sufficient to prove that Rev n {r) = 
in an interval < r < r G (n), where r (n) = 0(exp( zl\I^A ^y From Eq. (12.19) 
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one can prove that 

/ [f (r,x)] 2 V 1 (x)dx = Kr + 0(r 2 ), (12.29) 

Jo 

where K is a constant, K = 0(1). The integral does not vanish if r > 0. 

This suggests trying a double expansion in powers of g and r. However, again this 
will not lead to any restriction on Img n . The main problem is the relevant domain 
in t is small, i.e. r = 0(e~), and terms of order g 2 are much larger than terms of 
order r. 

To proceed further along the lines suggested by this paper one has to do two 
things: 

I. First, find an even function h(v), analytic for iseS(T ), and having no zeros in 
S(T ), such that if we define, x{ v ) as 

x (v) = £(v + l/2)h(v), (12.30) 

we have 

x (v)=0([t 2 ]-v), Kp<3/2. (12.31) 

where v = uo + it, t > T a . The point here is that x{ v ) is small but not smaller than 
0{g 2 ). 
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This first step is achievable. For example, we can define x as 
This will give 

X (i/) = 0([t 2 ]- (1+5) )- (12.33) 

The second requirement is much harder to achieve: 

II. One has to construct Jost functions, M^\iy, k), preferably of the Martin type, 
such that 

lim k ^oMM=x(v + l/2). (12.34) 

and for small g, ueS(T ), 

In addition has to be of the Martin type and it must have an asymptotic 
expansion in powers of g — [y 1 — 
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Appendix 

In this Appendix we first give a proof of the Laplace transform representation for 
the Marchenko kernel, F(i/;x), x > 0. 
Starting with the definition (8.2) 

F(v, x) = i- jf (5(1/; fc) - l^dfc, (A.l) 

where L is the line Imk = 5, with \ < 5 < 1, we note that for z/ e «S(T G ) we have 
^(z/; k) analytic in k for Imk > 5 > 0{e~^ 2 ~) ) except for the cut along the positive 
imaginary k-axis, 1 < Imk < oo. Second, in this region, we have a bound for large 
\k\ 

\S(v,k)-l\< ^,|fc|-oo. (A.2) 

This bound holds along any radial direction that excludes the cut. 

We can deform the contour L from along the line Imk = 1/4, to a contour 
surrounding the cut, i.e. 

F{y- x) = -!- f [S(is; k) - l]e tkx dk, x > 0. (A.3) 

where C starts at (— e, +ioo) and descends to (—e,+i), turns around the point k = 
i, and then extends from (+e,+i) to (+e,+ioo). The contribution from the large 
semicircle, \k\ = K, vanishes as K — > oo. For IK^ 1 ^ 2 < argk < n — \K\^ 1 / 2 ^ the 
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contribution is 0(exp(—\K\ 1 ^ 2 x)), and vanishes for x > 0. Here we use the bound 
(A. 2). For the regions, < argk < \K\~ 1 / 2 , and n — \K\~ 1 / 2 < argk < ir, the 
contribution from the semicircle to (A.l) will be 0(|i^| -1 / 2 ) and also vanishes as 
| If | — > oo. 

From (A. 3) we finally obtain 



where D is given in Eq. (8.5), and from Eq. (4.2) is just the discontinuity of S(u; k) 
along the cut. Noting that D = 0(e~ 7va ) as a — > oo, we see that (A. 4) will also hold 
for x = 0. 

The next task for this Appendix is to give a direct proof of the fact that U (u; x) 
and f^(y; k; x) as defined in Eqs. (8.93) and (8.94) do indeed satisfy the Schrodinger 
equation, i.e. to prove lemma (8.5) directly. We will also give an explicit expression 
for U(u; x). 

Following ref. 22 we define an operator, Q{y\ x), depending on two parameters, v 
and x, with Rex > 0, and v G S(T a ). Q acts on functions u(/3), 1 < (3 < oo, with 
u E L 2 (l, oo). We define Q as: 




(A.4) 




(A.5) 
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where 

1 Div 8)e~ 2f3x 

Q(u; x; a, (3) = , Rex > 0, (A.6) 

7r [a + p\ 

with D(v;(3) given by Eqs. (4.3) and (4.4). Q will have a finite Hilbert-Schmidt 
norm 

\Q(u;x)\ < Ae~ 2x , (A.7) 

where A depends on v. 

We introduce a new integral equation, 

Q(v;x;a,P)W(v;x;P)dp (A.8) 

with W e L 2 (l, oo). 

This integral equation is equivalent to the Marchenko equation (8.10). To see that 
we write 

Z(u; x; a) = -D(v; a)e- ax W(u; x; a) (A.9) 

A(v;x,y) is 

/oo 

Z(u;x;a)e~ ay da. (A.10) 

The Laplace transform exists since D = 0(e~ na ) and W G L 2 . Substituting (A.9) 
and (A. 10) in Eq. (A.8) and using (A.4) we obtain for A 

roo _ 

A(v; x, y) = F(is; x + y) + / A(v;x,u)F(v;u + y)a (A.ll) 
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But this is just the Marchenko equation which we have shown in section VIII does 
have a unique solution A. Hence A = A for all v e S(T ). Thus we conclude that Eq. 
(A. 8), which is of the Fredholm type, has a unique solution, since the homogeneous 
equation W = QW cannot have a solution for that will lead to the existence of a 
solution for A = AF which we have shown in section VIII is not possible. 

Given the function W{y \ x\ a) we can easily get expressions for U [y \ x) and / ± (^; k; x) 

The potential is given by 

U{u; x) = -2— [- / D(u, a)e' 2ax W(u; x; a)da],Rex > 0. (A.12) 
ax n Ji 

Similarly, from (A. 9) and (A. 10) we get 

k,x) = + e^(i f ^(^^)^^^^) ^, (A . 13) 

7r Ji (a ± ik) 

To check that we recover the same Jost functions we started with, we write 

7T Ji (ct ± «/eJ 

However the integral equation for H^ for a; = is trivially soluble. From Eq. (A. 8) 
we have 

ffM;a) = , + l/ %M^ ,. (A , 5) 

7r J (a + p) 

Setting 

0; a) = Af(->(i/; ia), (A.16) 
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and using the Eqs. (4.3) and (4.4) for D(u; a) we have 



M<-v; ia) = i + ^ - 1/4, r dp ma±m . (A , 7) 

Ji (p + a) 

This is our original expression for M^~\ 

The operator Q is a Fredholm type operator, and for Rex > 0, and v e «S(T ), we 
proved that there are no non-trivial solutions of the homogeneous equation u = Qu. 
Hence the determinant, Det(l — Q), cannot vanish for any Rex > 0. This determinant 
can be calculated explicitly, as was done in ref. 22, and we obtain 

Det (i-Q(,;,)) = i + i; R / ),n.kr^ 

(A.18) 

This series is absolutely convergent for all i?ex > 0, v e «S(T D ), and |imz/| < oo. 
Also as shown in ref. 22 the potential [/ (i/; x) is now given by Dyson's 25 formula, i.e. 

U(u;x) = -2^{log[Det(l - Q(u;x))]}. (A.19) 

Finally, we give here a direct check on the validity of the Schrodinger equation for 
^ and U (v\ x). 

From Eq. (A. 7) giving a bound on the norm of Q we see that for some x Q , \Q\ < 1 
for all (Rex) > x , x Q > Thus the iterative series for W(u;x;a) is absolutely 
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convergent for all Rex > x Q , 

oo 

w = 1 + Q n ( A -20) 

n=l 

Using Eqs. (A.6), (A.12), and (A.13), we obtain 

n=0 T 7l Jl (n j=0 («i + «i+l)) i=0 

(A.21) 

A similar series holds for f ± 

(A.22) 

Using these series we can check directly for Rex > x Q , that U and / ± give a 
potential and its unique Jost solutions. 
We define hP\ 

h (±) = e ±ikxj{±) (A>23) 

The Schrodinger equation for is now 
d 2 dh^ 

— ■/i ± (z/; fc; x) =F 2ik——(u; k;x) = U (u; x)/i ± (z/; fc; x) (A.24) 
ax 2 ax 

Substituting, the expressions (A.21) and (A.22) in the above, we see, after some 

algebra, that for x > x Q) (A.24) is satisfied if the following algebraic identity holds 

r— 1 n n r 

E{K + «n+l)(E %)}] = CC - a rCE a i)- ( A - 25 ) 

n=o j=o j=o j=o 
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But this is equivalent to 

r n r— 1 n r 

E «»(E«i) + EK+i)(E«;) = E«il 2 -- ( A - 26 ) 

ra=0 j=0 n=0 j=0 j=0 

This last equation is an identity and can be proved by induction. 

The Schrodinger equation thus is valid for all x > x Q . But again using analytic 
continuation, now in x we easily see that it must hold for all x > 0. All the terms in 
(A. 24) are analytic in the half plane Rex > 0. 
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Tables 



c(°) 


j = 










/ = 1 


1 










C(i) 


3 = 


1 2 








/ = 1 


6 


14 4 








2 


-14 


-8 








3 


8 










C(2) 


j = 


1 


2 


3 


4 


Z = 1 


36 


364 


500 


176 


16 


2 


-364 


-1000 


-528 


-64 




3 


1000 


1056 


192 






4 


-1056 


-384 








5 


384 











CP) 


j = 


1 


2 


3 


4 


5 


6 


/ = 1 


216 


7784 


29152 


29128 


10448 


1440 


64 


2 


-7784 


-58304 


-87384 


-41792 


-7200 


-384 




3 


58304 


174768 


125376 


28800 


1920 






4 


-174768 


-250752 


-86400 


-7680 








5 


250752 


172800 


23040 










6 


-172800 


-46080 












7 


46080 















Table 1: Values of C {n \l;j) for n < 3. 
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(n) 

Xi 


n= 1 


n = 2 


n = 3 


I = 1 


10.9973 


-460.8231 


28967.9828 


2 


4.7050 


-309.0434 


36560.3049 


3 


-7.4045 


451.5522 


-18626.4002 


4 


-8.2977 


910.6755 


-114291.8885 


5 




71.4580 


-76110.8714 


6 




-663.8194 


131495.9200 


7 






123526.6033 


8 






-111521.6508 



Table 2: Values of xT for n = 1, 2, 3. 

Figures 
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012345678 

Figure 1: Plot of V^x). 
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